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Abstract 

We discuss an extension of the (massless) Thirring model describing interacting 
fermions in one dimension which are coupled to phonons and where all interactions 
are local. This fermion-phonon model can be solved exactly by bosonization. We 
present a construction and solution of this model which is mathematically rigorous 
by treating it as a continuum limit of a Luttinger-phonon model. A self-contained 
account of the mathematical results underlying bosonization is included, together 
with complete proofs. 


1 Introduction 

Exactly solvable models have played an important role in the history of condensed mat¬ 
ter physics. A famous example is the Luttinger model [24, 37, 38] solved in an important 
paper by Mattis and Lieb [27]. The key to this solution is a collection of mathematical 
results commonly known as bosonization (another common name is boson-fermion cor¬ 
respondence). In the case of the Luttinger model, bosonization is particularly powerful: 
not only all energy eigenstates and eigenvalues but any quantity of physical interest can 
be computed exactly by analytical means; see [18] and references therein. Moreover, as 
shown in [25], a non-trivial limit of this solution of the Luttinger model is equivalent to 
Johnson’s solution [19] of the massless Thirring model [37]. 

Bosonization is today a standard tool in condensed matter physics, and it has allowed 
physicists to obtain interesting results even in situations where exact solutions are not 
available. However, it is not always easy to judge the mathematical status of such results 
in the physics literature. One difficulty is that bosonization has a complicated history, 
and papers developing its mathematical foundations were written with a focus on other 
applications, notably the representation theory of infinite dimensional Lie algebras and 
string theory; see e.g. [4, 6, 13, 20, 29]. There exist excellent reviews and textbooks on 
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bosonization in the condensed matter physics literature, including [11, 14, 15, 17, 22, 32, 
33, 35, 36, 39, 41], which provide the reader with practical knowledge on how bosonization 
is used in applications. However, to our knowledge, the existing literature lacks a self- 
contained account of all mathematical details needed to promote formal bosonization 
results in condensed matter physics to rigorous mathematical statements. One aim of 
the present paper is to fill this gap, and for this we combine different approaches found 
in the mathematics literature [6, 20, 34]. 

We discuss a model of interacting fermions coupled to acoustic phonons which can be 
solved exactly by bosonization. Our program is to develop all mathematical tools which 
are needed to give this result, known in the condensed matter physics literature for a 
long time [8, 12, 40], a rigorous mathematical foundation. In particular, we include a 
self-contained account of bosonization, where we try to remain close to how bosonization 
is used in condensed matter physics, but, at the same time, also address the relevant 
functional analytic issues. 

The fermion-phonon model we consider is (formally) defined on the real line M by 
the Hamiltonian^ 

H= f dx(vF'^:ipl{x)r{-idx)iJr{x):+^:(llp{xf + Vp[dx^p{x)fy. 

+ X:ip\_{x)'ip_^_{x): :ip^_{x)'ip_{x): + g E :'il)l{x)i;^{x):da;^p{x)j (1.1) 

r=± ^ 

with fermion fields satisfying canonical anticommutation relations (CAR): 

= ^r,r'S{x - y), ^= 0’ (1-2) 

and phonons described by Hermitian boson fields 4>p(x) and np(x) satisfying canonical 
commutation relations (CCR): 

[^p{x),Up{y)] = i5{x - y), [4>p(x), 4>p(y)] = [np(x), np(y)] = 0, (1.3) 

and which commute with the fermion fields. The colons in (1.1) indicate (Wick) normal 
ordering. The fermion-phonon model has four parameters: the Fermi velocity vp > 0, 
the phonon velocity up > 0 with vp < vp (typically vp <C vp), the fermion-fermion 
coupling constant A, and the fermion-phonon coupling constant g. As will be shown, if 
the coupling constants are in the range 

X<2pvp, 2 {g/vp)"^ < 2pvp + X, (1-4) 

then this model can be constructed and solved by bosonization. (If these conditions are 
not fulfilled, the model describes an unstable system.) 

We note that, in our approach to the fermion-phonon model, we follow previous work 
on the construction and exact solution of a two-dimensional interacting fermion model 

^The notation used here is common in the physics literature but suppresses mathematical details. 
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by bosonization in [9]; see [10] for a review. However, these papers only stated the 
mathematical results underlying bosonization, and thus one motivation for the present 
paper is to provide self-contained proofs of these results. 

The fermion-phonon model in (!.!)-(1.4) is interesting since it has zero-temperature 
phase transitions to charge-density-wave and superconducting states, and its phase di¬ 
agram can be determined from certain fermion four-point correlation functions [8, 40]. 
From a mathematical point of view, this model is an extension of the massless Thirring 
model, which (essentially) corresponds to the limiting case g = 0- However, as opposed 
to this limiting case, the fermion-phonon model is not relativistically invariant, which is 
obvious since its correlation functions depend on two different velocities vp and vp] see 
(1.5) below. Nevertheless, as will be discussed, this model shares many properties with 
relativistically invariant quantum held theory models. Since many effective low-energy 
models in condensed matter physics happen to be relativistically invariant, this is an 
interesting counterexample. 

As a further motivation of the bosonization techniques presented in this paper, we 
mention recent mathematical work on bosonization based on functional integral tech¬ 
niques; see e.g. [2] and references therein. These methods allow for treating even models 
which need not be exactly solvable, but they have been restricted, up to now, to models 
with fermions only. It would be interesting to extend these results to models where 
fermions are coupled to bosons. We hope that the results in the present paper can serve 
as a useful guide in this direction. 

As will be shown, the fermion-phonon model can be made mathematically precise 
by introducing cutoffs L > 0 and a > 0, which have the physical interpretation as sys¬ 
tem size and interaction range, respectively. The regularized model thus obtained is an 
exactly solvable extension of the Luttinger model. To obtain the correlation functions 
for the model in (1.1)-(1.4), one has to pass to both the continuum limit a —)■ O"'' and 
the thermodynamic limit L —)■ oo; this is a generalization of the limit in [25] where the 
Thirring model is constructed from the Luttinger model. As will be shown, the contin¬ 
uum limit is non-trivial in that it requires additive and multiplicative renormalizations. 
However, passing to these limits simplifies the correlation functions significantly, and 
physical arguments suggest that the result is universal in the following sense: the phys¬ 
ical properties of the fermion-phonon model at intermediate distances do not depend 
on the regularization, i.e. on short-distance details. For example, after multiplicative 
renormalization, the fermion two-point correlation functions we obtain in these limits 
are 


(V'^(t,x)V’]t/(0,0)) = 




n 

X=F,P 


U 

rx — vxt -|- i0+ _ 


\Px 


it 

—rx — Vxt + 


4 


(1.5) 


with parameters vp, ap and vp, pp, Up determined by the model parameters (see 
Proposition 4.4); > 0 is an arbitrary length parameter introduced by the multiplicative 

renormalization. 

The rest of the paper is organized as follows. The collection of mathematical results 
known as bosonization are presented and proved in Sections 2 and 3. For this we use an 
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infrared cutoff L > 0 corresponding to the length of space, i.e. the real line M is replaced 
by the interval [—L/2,L/2]. These results are used in Section 4 to construct and solve 
the fermion-phonon model described above. Section 5 contains concluding remarks. In 
order to abbreviate our formulas certain variables are reserved to belong to specific sets 
as given by the following table: 


Variables 

Set 

r, r' 

{+) “} 

X, x' 

[-L/2,L/2] 

k,k' 

(27r/L)(Z + l/2) 

p,p' 

(27r/L)Z 

X, X' 

{F.P} 


We also use the convention that sums over variables in this table range over all allowed 
values, unless specified otherwise. These and other conventions are summarized in Ap¬ 
pendix A. In Appendix B we give a summary of the results in the paper using a notation 
common in the physics literature, in order to connect our results with the latter. It 
also contains comments on the history of bosonization. Results from functional analysis 
which we need are collected in Appendix C. 

2 From fermions to bosons 

The results referred to as bosonization state that two different quantum field theory 
models, one describing fermions and the other bosons, are actually the same. There are 
two approaches to this: one can start with the fermions and construct the bosons, or one 
can start with the bosons and construct the fermions (see Figure 1). Regardless of this 
choice, the same collection of results is obtained in the end. In our opinion, the second 
approach is simpler from a mathematical point of view. The first approach is, however, 
more natural from the point of view of condensed matter physics, and this is why we 
follow it in this paper; see e.g. [5, 35] for the second approach. 



Figure 1: Symbolic illustration of bosonization. 

The present section is concerned with the construction of the bosons from the 
fermions; the reconstruction of the fermions from the bosons is addressed in Section 3. 
We remind the reader that some notational conventions used throughout this paper are 
explained in Appendix A. 
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2.1 Fermion Fock space and Hamiltonian 

Our first aim is to construct a Hilbert space T^ the so-called fermion Fock space, giv¬ 
ing a precise mathematical meaning to the quantum field theory model defined by the 
Hamiltonian^ 

(2.1) 

r k 

for r = ± and k G {27rjL)^!^ 1/2), with fermion field operators iprik) and = 

ifrik)^ characterized by the CAR 

= ^Sr,r'6k,k', (^ 0 } = 0 - ( 2 - 2 ) 

The physical requirement underlying this construction is that Hq must have a ground 
state. This is ensured by requiring that there is a distinguished element H in F, the 
so-called vacuum, which satisfies 

= 'tpl{—rk)^ = 0 VA: > 0 (2.3) 

and is normalized, i.e. (H, H) = 1, where (•, •) is the inner product in F. The Hilbert space 
F is fully determined by these conditions, and it is such that FLq defines a self-adjoint 
operator on F bounded from below. We now recall how this can be substantiated. 

To construct F we use fermion creation operators cl{k) satisfying the conditions; 

|c^(A:),c|;,(/c')| = 5r,r'4,fc', = 0, c^(/c)H = 0, (2.4) 

with Cr{k) = cl{k)^ and a normalized state Q ^ F. Assuming that ci^\k) are operators 
on F, one can construct other states in F as follows: 

c\{ki) ■ ■ ■ c^^{kNjf)(^-{k'i) ■ ■ ■ c^_{k'j^_)^ for /ci > ... >/cAr+j k'i> ... > k'jq_, (2.5) 

where N± are non-negative integers. Using the relations in (2.4), it is straightforward to 
prove by induction over A^_|_ and N- that these states are orthonormal. Let Vp denote 
the vector space of all finite linear combinations of the states in (2.5). This set forms a 
pre-Hilbert space, and F is obtained from T>f by norm-completion; see e.g. [30]. 

We find it convenient to introduce the following notation for the states in (2.5): 

Vn = (nct(A:)-’'W JH, nrik) G {0,1}, (2.6) 

\r,k / 

where only finitely many of the quantum numbers nr{k) are non-zero and n denotes the 
infinite vector {nr{k))^^. The arrow over the product symbol indicates that the fermion 
operators are ordered as in (2.5). 

^This Hamiltonian corresponds to an important model in particle physics describing relativistic 
fermions in 1-1-1 dimensions; cf. Remark B.l in Appendix B. 
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The above definitions allow us to write formulas in a compact way. For instance, 
orthonormality and completeness of the states in (2.5) can be expressed as 

{Vn,Vn') = Vn iVn , ■) = I ■ (2-7) 

n 


Note that the states are common eigenstates of cl{k)cr{k), i.e. 

4ik)Cr{k)Vn = nr{k)r]n ( 2 - 8 ) 

(this follows from the CAR in (2.4) and the definition in (2.6)). 

Remark 2.1. It is not difficult to construct a bijection between the set of quantum 
numbers n and N to prove that the Hilbert space J- thus constructed is separable. 0 

We are now ready to give a precise meaning to the model defined in (2.1)-(2.3). The 
method we employ is construction by exact solution in the sense that self-adjointness of 
the Hamiltonian is shown by constructing all its eigenstates and eigenvalues. 

Using the relations in (2.4) it is easy to check that the operators 


-Iprik) = 


[17 (i.\ 


if rk > 0, 
if rk < 0, 


YJ(/c) = Yr(^)^ 


(2.9) 


satisfy the relations in (2.2)-(2.3). Moreover, for the operators we consider, normal 
ordering can be made precise as follows: 


lA: =A-{n,An). 


( 2 . 10 ) 


This implies 

:ipl{k)[r{k)- = sgn{rk)^cl{k)c^{k), ( 2 . 11 ) 

ZTT 

which allows us to cast the Hamiltonian in (2.1) in a manifestly positive-definite form: 

Ho = Y^ '^\k\cl{k)c^{k). ( 2 . 12 ) 

r k 

Remark 2.2. The CAR in (2.4) and the eigenvalue equation (2.8) imply that the fermion 
operators c^\k) are bounded on F with operator norm equal to unity, and thus cj(/c) 
is the Hilbert space adjoint of Cr{k). It follows that [[4\k) are also bounded operators. 
This implies that the algebra generated by the fermion field algebra is a (7*-algebra, 
which is the starting point of an elegant mathematical approach to fermion quantum 
field theories; see e.g. [3]. However, to remain self-contained, we follow a more pedestrian 
approach avoiding C^-algebra techniques. 0 

From the equations above one can deduce several important properties of Hq, which 

we summarize as follows. 
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Lemma 2.3. The Hamiltonian Hq in (2.1) defines a self-adjoint operator on T, and it 
satisfies the following commutation relation onVp: 

HQ,'4>l{k) =rkfil{k). (2-13) 

Moreover, all states in (2.6) are eigenstates of Hq with corresponding eigenvalues 

(2.14) 

r k 

In particular, Q is the ground state of Hq with ground state energy zero: HqQ = 0. 


Proof of Lemma 2.3. The eigenvalue equation 


77o??n = -^n^n 


(2.15) 


follows from (2.8) and (2.12). This implies that Hq is a self-adjoint operator (this is a 
consequence of a simple general result stated in Lemma C.3). 

We note that the commutation relation in (2.13) is easily motivated by formal com¬ 
putations using the definition of Hq in (2.1) and the CAR in (2.2) if we ignore normal 
ordering. To make these computations precise we introduce a cutoff A > 0 and study 
the regularized Hamiltonian 


:fil{h)'ipfik): 0(A - |/c|) 

r k 

Ott 

= ^ ^ -—rki)l{k)fifik) + Ea, 


(2.16) 


where Ea = (L/(27r))A^ -|- 0(A) is a positive constant which diverges with increasing A. 
The second equality is obtained using ■.fil{k)fir{k): = ^J(/c)'0r(^) — (T/(27r))0(— r/c), 
where 0(-) is the Heaviside function (this is a simple consequence of the definitions). 

Since the regularized Hamiltonian Hq is a finite sum of bounded operators, it is 
bounded, and thus commutation relations, etc., can be safely derived by formal compu¬ 
tations using the CAR in (2.2). In particular, it is obvious that limA-s-oo 77,)^ = Hq on 
Dp (in a sense made precise in Definition C.5). 

We now observe that, for each state in (2.6), 


HQ,'4>l{k) r/n = HQ,'4>l{k) 


(2.17) 


independent of A for sufficiently large values. It follows that the result obtained by 
formal computations is correct if interpreted as an identity on the domain T>p (this 
follows from a simple general result stated in Lemma C.6). □ 
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Remark 2.4. Throughout this paper we use the same strategy to establish commutation 
relations for unbounded operators A: we introduce a cutoff A > 0 and approximate A 
by a family of bounded operators {A^}a>o such that A = limA-s-oo A^. This strategy 
is simple in that it avoids topological arguments to make precise mathematical sense of 
such limits. To not have to repeat the same technicalities each time we use this strategy, 
we state the precise definition of this limit, together with a simple abstract lemma, in 
Appendix C; see Definition C.5 and Lemma C.6. 0 


2.2 Fermion densities 

We now introduce the fermion densities 

Mp) ■^Iik)'iprik + P)- (2.18) 

k 

for p G (27r/L)Z, which we regard as a definition motivated by the Fourier transform of 
the (formal) fermion densities in the introduction (cf. (B.3) and (B.12)). 

Proposition 2.5. The fermion densities Jr{p) in (2.18) are well-defined operators on 
T mapping Vp to Vp. Moreover, they satisfy 

JAp)^ = Jri-p)^ (2-19) 

X{rp)n = 0 Vp > 0, (2.20) 


together with the following commutation relations onVp: 


Hq, Jr{p) 
Jr{p),Jr'{p') 


6ryAl{k-p), 

-rpJr{p), 

— — Op^—pf . 

^71 


(2.21a) 

(2.21b) 

(2.21c) 


In particular, the zero modes Triff) are self-adjoint operators, and all states in (2.6) are 
eigenstates of Jr (0) with corresponding eigenvalues 


qr^n ^ sgn(rfc)n)-(fc). (2.22) 

k 

Proof. Using (2.9) one can write 

MO) = E sgn(r/c)cj.(/c)c^(/c), (2.23) 

k 


and thus (2.8) imply that the states in (2.6) are eigenstates of Jr(0) with eigenvalues 
given by (2.22). This also proves that Jr(0) are self-adjoint operators (cf. Lemma C.3). 
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It is easy to motivate (2.19) and the commutation relations in (2.21a)-(2.21b) by 
formal computations using and the CAR in (2.2). However, the commu¬ 

tation relation in (2.21c) is non-trivial (a formal computation suggests that it is zero), 
and this highlights the need for precise arguments. 

A potential difficulty with the fermion densities is that the sum on the r.h.s. in (2.18) 
is infinite. As for Hq in Section 2.1 (recall Remark 2.4), we deal with this by introducing 
regularized fermion densities 

^ _ O-jr 

JHp) = + p): 0(A- |A:-b II), (2.24) 

k 

where A > 0 is a cutoff and ©(•) is the Heaviside function. Note that these operators 
are bounded for A < oo by the same argument as for Hq in the proof of Lemma 2.3. 
We use (2.9) to write, for all p > 0, 

X{rp) = c\.{k)c^{k + rp) + ^ c^{k)c^{k + rp) — ^ c\{k + rp)c^{k)^ 

rk>0 —p<r/c<0 rk<—p 

Jr{—rp) = ^ c\{k)c^{k — rp) + ^ c\{k)c\{k — rp) — ^ c\{k — rp)c^{k). 

rk>p 0<rk<p rk<0 

(2.25) 

This and (2.23) make manifest that all Jrip) are well-defined operators mapping Dp to 
Vp. It also proves (2.20) and implies that, for each state t]^ in (2.6), 

Jr{p)Vn = JriP)Vn (2-26) 

independent of A for sufficiently large values. Thus, limA-s-oo Jr{p) = Jr{p) on Vp (in 
the sense of Definition C.5). Since it follows from the definitions that J^{p)^ = J^{—p), 
this implies (2.19). 

Straightforward computations using the CAR in (2.2) imply 

Jr{p),^l{k) =Sr,r'ipl{k-p)QiA-\k-^\). (2.27) 

Acting with this identity on a fixed state rj^, for fixed k and p, the r.h.s. is independent 
of A for sufficiently large values. This proves the commutation relation in (2.21a) on T>p 
(cf. Lemma C.6). 

To prove (2.21b) we use the CAR in (2.2) to compute 

Ho,Jr'{p) ='^‘^r{kA-{k + p)A)^l{k)'4}^{k + p)e{A'-\k + l\), (2.28) 

k 

where /ca = kQ{A — |/c|). As a simple consequence of the definition of normal ordering, 

i’lik)^rik') = ■'4’lik)i’rik')-+ ^h,k'Q{-rk). 


(2.29) 
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It follows that 'tpl{k)il)r{k+p) on the r.h.s. of (2.28) can be replaced by :'ipl{k)ipr{k + p)- 
(this is obvious for p ^ 0, and it is also true for p = 0 since then the r.h.s. is identically 
zero anyway). Acting with the identity thus obtained on a fixed state for fixed p, 
the r.h.s. becomes independent of A and A' when both are sufficiently large, and this 
implies (2.21b) (cf. Lemma C.6). 

We are left to prove (2.21c). The case r ^ r' is trivial, and we therefore only consider 
the case r = r'. By straightforward computations. 


Jr{p), Jr (P) — X/ l^+2l)®(^^“ l^+P+^l) 


k 

- 0(A' -\k + |-|)0(A - \k+p + ^\)^ipl{k)'il>j.{k + p + p). 


(2.30) 


The r.h.s. can be cast in a normal-ordered form by inserting (2.29). It is thus equal to a 
sum of two terms; the first is equal to the r.h.s. in (2.30) but with il>l{k)'4>r{k + p + p') 
replaced by :il)l{k)'tpr{k + p + p'): , and the second is the C-number term coming from 
normal ordering. Acting with the identity thus obtained on a fixed state we find 
that, for fixed p and p' and sufficiently large A and A', the contribution of the first term 
is independent of the cutoffs and thus zero. The only remaining contribution is therefore 
the second term; 

JHp), Jr'{P) 

k 


Vn ^p+p',0 0(A - + II) - 0(A - - II) 0(-rfe)ry^ (2.31) 


assuming A' > A, without loss of generality. For sufficiently large A, the r.h.s. in (2.31) 
is independent of A and equal to 5pj^p'^Q{rpL/{27r))ri^. This implies the result (cf. 
Lemma C.6). □ 


2.3 Bosons, charges, and Klein factors 

In this section, the physical interpretation of the fermion densities Jr{p) as bosons (for 
p 7 ^ 0) and chiral charge operators (for p = 0) is explained. We also discuss how this 
motivates the introduction of the Klein factors R± as charge-changing operators. 
Consider operators b^^\p) defined for p / 0 as follows; 


b{p) = 



if p > 0, 
if p < 0, 


b^p) = 6(p)^ 


(2.32) 


Proposition 2.5 implies that these are well-defined and obey standard relations for boson 
annihilation and creation operators on 


b{p),b^{p) 


6. 


p,p'j 


b{p),b{p) 


= 0 , 


6(p)n = 0. 


(2.33) 


The zero modes J± (0) play a special role as chiral charge operators, and the notation 

Q± = J±(0) (2.34) 
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will sometimes be used to emphasize this. Indeed, they are self-adjoint operators with in¬ 
teger eigenvalues (cf. Proposition 2.5) which, as discussed below, have a natural physical 
interpretation as chiral charges. 

Our aim is to generate the fermion Fock space J- from the vacuum using bosons. 
However, since the zero modes Q± commute with the boson operators b^'^\p) and since 
= 0 (cf. (2.21c) and (2.20)), the boson operators in (2.32) are not enough: they 
can only generate states t] with charge zero {Q±r] = 0). This problem is solved by the 
Klein factors R ±, which are unitary operators on T commuting with all boson operators 
but not with the zero modes: 


QrRr' = Rr'iQr + (2.35) 

This implies that, if r/ is a zero-charge state {Q±ri = 0), then R_^~r} is a state with 
charges (g+, ?-): 

Q^R'‘+RZ^^-ri = q^R^^^ RZ"^-T] (2.36) 

for all integer pairs (g+, g_), where RZ'^ = integers q. 

Remark 2.6. The picture above can be summarized as follows: 

R= 0 (2.37) 

(g+,g_)eZ2 

where is the subspace of R containing all states rj such that Qrrj = qrf], and the 

unitary operator RZ^Rj^~ maps the zero-charge subspace bijectively to 

The subspaces are known as charge sectors or superselection sectors. 0 

A construction of the Klein factors, together with proofs of properties needed in 
practical computations, is given in Section 2.4; see Proposition 2.7. One important 
result underlying bosonization is that the boson operators {p) , together with the Klein 
factors R±, generate R] see Proposition 2.12. 

In the remainder of this section we elaborate on the physical interpretation of Q± as 
chiral charge operators. This interpretation is based on Dirac’s hole theory and can be 
skipped without loss of continuity. 

The model in (2.1)-(2.3) describes a system of fermions with one-particle states 
labeled by pairs of quantum numbers (r, k) and with one-particle energies rk. The many- 
body ground state D is a Dirac sea where all one-particle states with negative energy 
are filled and all one-particle states with positive energy are empty; this is expressed by 
(2.3), where the second relation is a consequence of the Pauli principle. The physical 
interpretation of (2.9) is that Ylik) is a creation operator for a particle if rk > 0 and an 
annihilation operator for a hole if rk < 0. Similarly, V’r(^) is an annihilation operator 
for a particle if rfc > 0 and a creation operator for a hole if rk < 0. Moreover, the 
quantum numbers nr{k) used to label the states in (2.6) correspond to the excitation 
numbers of the one-particle states {r,k); the excitations are particles for rfc > 0 and 
holes for rk < 0. Finally, cl.{k)cr{k) are excitation-number operators with eigenvalues 
equal to either 0 or 1 (this corresponds to the Pauli principle). It is natural, from a 
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physical point of view, to assign charge numbers +1 to particles and —1 to holes (Dirac 
proposed this interpretation in a particle physics context where the holes correspond to 
antiparticles). This suggests to interpret Qr as an operator which measures the charge 
of particles and holes with chirality r (cf. (2.23) and (2.34)). 

It is also worth to mention a beautiful mathematical interpretation of the fermion 
densities and the Klein factors as implementers of gauge transformations. This inter¬ 
pretation clarifies that the Klein factors exist due to topological reasons; see the final 
paragraph in Section 3.1. 

2.4 Construction of the Klein factors 

This section contains a construction of the Klein factors R±. We start with a summary 
of the main results: 

Proposition 2.7. (Klein factors) There exist unitary operators R± on J- mapping 
Vp to Vp and satisfying the eommutation relations 

Jr{p),Rr' = rdrySpflRr, (2.38) 

Ho,Rr = r^^Jr{0),Rr'^ (2.39) 

on Vp, together with 

R+R- = -R-R+. (2.40) 

(A constructive proof is given further below.) 

Note that (2.38) summarizes properties of the Klein factors mentioned in Section 2.3; 
choosing p ^ 0, it implies that R± commute with the boson operators in (2.32), and 
choosing p = 0, it implies (2.35) (recall that Qr = Jr(0)). 

The following corollary of Proposition 2.7 contains a list of identities which we state 
for later reference; 

Corollary 2.8. For all q± £ h and c G C, the following identities hold true: 


R^^+Rj- = Rj-R‘^+, (2.41a) 

Qr,RfRZ'^-]=qrR'^^RZ‘^-, (2.41b) 

^lcQrR±q± ^ ^tCqrSr,± ^IcQr (2.41c) 

on Vp, together with: 

QrR'^+RZ‘^~^ = qrR'^+RZ‘^~^, (2.42a) 

HoR'^_^ZRZ‘^-n = y(g+ + q^_)R'^+RZ‘^-n, (2.42b) 

IJ 


= <5,^,o<5g_,o. 


and 


(2.42c) 
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(The proofs of these identities are straightforward and thus omitted.) 

It is interesting to note that the Klein factors can be fully characterized by the 
following conditions: 

Rj^2\k) = + ^)R±, = -i’2\k)R±, (2.43) 

R^n = R-^ = ^_(f )fl. (2.44) 

In our proof of Proposition 2.7 below we show that these conditions define unique oper¬ 
ators with the stated properties. 

Remark 2.9. Equation (2.43) is a simple example of a defining condition for an im- 
plementer of a Bogoliubov transformation; see [31] for a unified construction of such 
implementers. 0 

Remark 2.10. The conditions in (2.43)~(2.44) are not minimal [31] and are chosen so as 
to allow for a simple proof. 0 

In our proof of Proposition 2.7 we make repeated use of the following technical result: 

Lemma 2.11. Let C be a linear operator from Vp to Vp satisfying the following two 
eonditions: C commutes with all fermion field operators 'il:2\k), and CQ = cQ for some 
eonstant c G C. Then C = cl. 

Proof. The first condition and (2.9) imply that C commutes with all operators c2\k). 
The action of C* on a fermion state in (2.6) is then fully determined by its action on 
n, and the second condition then implies that Ctj^ = crj^. The assertion follows from 
linearity. □ 

Proof of Proposition 2. 7. We rewrite the defining conditions of R± in terms of the op¬ 
erators c2\k) using (2.9). The identities in (2.43) take the following form: 

if= 
otherwise, 

R^c^flik) = -p2l{k)R^, (2.45) 

if= 
otherwise, 

= cl{l)n. (2.46) 

Using this, the action of R± on the fermion states in (2.6) can be computed, and the 
result is obviously in Dp. It then follows that (2.43)-(2.44) determine two unique linear 
operators R± mapping Vp to Vp. 

We now show that both RrRt and rIRt- satisfy the conditions of the operator C 
in Lemma 2.11 with the constant c = 1, which then implies unitarity: rI = Rf^ (cf. 
Lemma C.4). By taking the adjoint of (2.43) and replacing k with k — ^ we obtain 

Ri22\k) = - x)4> 444) = -44)4> 



_ J 

Cr{k+^)R, 
^r{j-)R 


cl{k+‘^)Rr 


(2.47) 
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which together with (2.43) implies that RtB}- and E^-Rr commute with all fermion field 
operators To check the second condition, we compute the inner product of 

with the fermion states in (2.6). Using (2.45), (2.46), and the CAR in (2.4) we find 


(r?^,i?tU) = (R,r?^,U) 


1 if r?^ = cJ(-f)U, 
0 otherwise 


(2.48) 


(to see this, we observe that = Cr{^)Rr^ = Cr(J)cJ(J)n = and that no 

other inner product can yield a non-zero result), which is only possible if 

Rln = cl{-l)Q. (2.49) 

It follows that rIRt^ = i?jcj(-^)n = Cr(—^) rJu = and similarly RrRl^} = fl. This 
concludes the proof of the unitarity of the Klein factors. 

We now turn to the commutation relation in (2.38). Using the defining properties of 
the Klein factors in (2.43) and the regularized fermion densities in (2.24) to compute 

(r,, J('(p)4 - J,^(p)) r? (2.50) 


for an arbitrary r] one finds that this is equal to 

^ry X] - 1^ - X + il) “ - 1^ + + p)r]. (2.51) 

k 

The fermion bilinear above can be cast in a normal-ordered form by inserting (2.29), 
and, by a similar argument as the one following (2.31), we conclude that (2.51) is equal 
to 

5ry ^(^0(A — ^ -I- ||) — 0(A — \k + ||)^0(—r/c)r/ = —rSryr] (2.52) 

k 

for sufficiently large A. This proves RpJ^{p)I^p — J^{p) = —rS^y on Vp (cf. Lemma C.6), 
which is equivalent to (2.38). 

The commutation relation in (2.39) can be proved in a similar manner by showing 
that Hq — Rrli^Rl- = r{p/L){Jr{R) + RrX(O)iiJ) on Vp. For this we use (2.16), (2.24), 
and (2.29) to compute 



r^h^rI 


TT 



J(^(0) + J^^(0)4 


T] 



1) (0(A -\k\)- 0(A -\k- ^f\))i^l{k)Mk) 


+ r^(e{A - \k - ^\) + e{A - \k\)^Q{-rk)\r] (2.53) 



2 From fermions to bosons 


15 


for an arbitrary r] ^T)p. As before, the fermion bilinear can be cast in a normal-ordered 
form by inserting (2.29), and, by a similar argument as the one following (2.31), we find 
that the r.h.s. in (2.53) is equal to 

r ^(A:0(A - |A:|) - {k - ^) 0(A - \k - ^\)^e{-rk)r] = 0 (2.54) 

k 

for sufficiently large A, which proves the result (cf. Lemma C.6). 

Finally, equation (2.40) can be proved by showing that satisfies the 

conditions of the operator C in Lemma 2.11 with the constant c = — 1. This implies 
that R_^_R_R}^R}_ = —I, and the result then follows by multiplication with R-Rp from 
the right. Indeed, it is easily verified using (2.43) and (2.47) that Rj^R_R}j^R^_ commutes 
with all fermion field operators 'ip^\k). Moreover, (2.45)-(2.46) and (2.49) imply that 
Rj^R_Ftj^R}_VL = c_( j)c_,_( J)cl( j)c^^_( J)n = —fl, where the CAR in (2.4) was used in 
the last step. □ 

2.5 Boson-fermion correspondence 

We are now ready to explain in which sense the fermion Fock space J- is identical with a 
boson Fock space. For this we recall the boson creation and annihilation operators b\p) 
and h{p) in (2.32). Since the Klein factors commute with the boson operators (this is a 
simple consequence of (2.38)), all states R^j^ where q± G Z, constitute a possible 
boson vacuum: 

6(p)R)(.+ R:®-fI = 0 Vp/0. (2.55) 

This suggests to introduce boson states 

= 1 n I R^R-J-^, q± G z, m{p) G Nq, (2.56) 

\p^o VMp)'- ) 

where m is short for the vector ((m(p))p^Q ,g+,g_) and where only a finite number of 
the quantum numbers m{p) are non-zero. 

Our main result in this section is that the boson states in (2.56) comprise a complete 
orthonormal basis in the fermion Fock space J-: 

^ vSiivL ■) = ^■ (2-57) 

m 

Orthonormality is easy to check (it follows from the commutation relations in (2.33), 
the vacuum condition in (2.55), and the unitarity of the Klein factors). The non-trivial 
part is therefore completeness. 

Recall that Dp the vector space of all finite linear combinations of the fermion 
states in (2.6). Since the boson creation operators {p) and the Klein factors map Rp 
to T>p, it is clear that all boson states in (2.56) are in T>p. The converse is also true: 
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Proposition 2.12. The vector space of all finite linear combinations of the boson states 
in (2.56) is identical with Vp. 

Our proof of Proposition 2.12 is based on the following: 

Lemma 2.13. All boson states rj^ in (2.56) are eigenstates of the fermion Hamiltonian 
Hq in (2.1) with corresponding eigenvalues 

E^^2{(ll + <l-) + ^\p\'m{p). (2.58) 

p^O 

Proof. The definition of the boson operators b^^\p) in (2.32) and the commutation re¬ 
lations in (2.21b) imply 


Ho,b^{p) 


\p\h^{p) Vp / 0 


(2.59) 


on Vp. This, together with (2.42b), yields the eigenvalue equation HqT]^ = 

Proof of Proposition 2.12. Let E he a. fixed eigenvalue of Hq. Denote by dimj7’(£') the 
number of distinct fermion eigenstates such that = E (cf. Lemma 2.3). Corre¬ 
spondingly, denote by dimB(i?) the number of distinct boson eigenstates i]^ such that 
E^ = E (cf. Lemma 2.13). It follows from (2.14) and (2.58) that both dimi?(i?) and 
dimB(Fi) are finite. 

We now prove that 

dimi7’(£') = dims(F^) (2.60) 

for all eigenvalues E of Hq. To this end, we introduce the generating functions 

Zj = ^ <Hmx{E)e-^^ {X = B, E) (2.61) 

E 

for some parameter /3 > 0, where the sum is over all distinct eigenvalues E of Hq. (As 
the symbols suggest, fi and have physical interpretations as inverse temperature and 
partition functions, respectively.) 

The generating functions can be computed as follows: 

4 = = (2.62) 
n m 

From the explicit expression for the eigenvalues in (2.14) we obtain 

r k nr(fc)e{0,l} 

nSZ / Vn=l 



4 


(2.63) 
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and similarly, using (2.58), 


^B= n 


q'_l_,g'_GZp/0 m(p)=0 
/ oo oo 


E 


n E ^ 

n=l m=0 


Hence, Zp = Zo is equivalent to 


fl27r \ 
—p^nm I _ 




(2.64) 


\q£Z 


o 277 




oo - 

2 \ / 1 —r 1 


\n=l 


E^’ )(n 


n=l 


1 - 


(2.65) 


with z = . The latter is a special case of a mathematical identity known as 

Jacobi’s triple product; see e.g. Chapter 16 in [1]. It follows that = Z^ holds true, 
which proves (2.60). 

In conclusion, we have proved that every fermion state can be written as a linear 
combination of boson states 77 ® such that = E^, and there is only a finite number 
of such boson states. The converse is obvious (see the above). □ 


2.6 Boson Hamiltonian 

From the results in Section 2.5 it is easy to deduce that the fermion Hamiltonian Hq in 
(2.1) can be written as a boson Hamiltonian: the commutation relations in (2.33) and 
(2.42b) imply that the boson states in (2.56) are eigenstates of the operator 

Hi = I (<5+ + Q-) + (2-66) 

with corresponding eigenvalues E^^ in (2.58). Thus, due to Lemma 2.13, it follows that 
{Hq — Hi)rii^ = 0 for all boson states. Since these states form a complete orthonormal 
basis in E, we conclude that Hq = Hq on E. Expressing the operators (p) and Q± in 
terms of the fermion densities (cf. (2.32) and (2.34)) we can state this result as follows. 

Proposition 2.14. (Kronig’s identity) The following holds true as an identity of 
self-adjoint operators on E: 

r \ p>0 / r p 

Remark 2.15. Mathematical physicists often write Kronig’s identity in (2.67) with boson 
normal ordering, i.e. with iJr{—p)Jr{p)i instead of :Jr{—p)Jr{p)-- We note that for 
expressions containing only factors with at most two field operators, as in this case, both 
normal orderings give the same result. However, for other expressions it is important to 
distinguish between boson and fermion normal ordering. 0 
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We conclude this section by making precise how this is related to a conventional 
boson Hamiltonian. It is natural to introduce the following operators; 


n(p) = (j+ip) - J-{p)), 

^(P) = [-Uip) + J-ip)) 

which obey CCR of neutral bosons: 







^{p),^Hp') 


{p / 0 ) 


flip),&{p') 


= 0 , 


( 2 . 68 ) 


(2.69) 


where n^(p) = n(—p) and $^(p) = ^{—p). Furthermore, they allow Kronig’s identity to 
be written in the form^ 

Ho = j {Ql + Q-) + E Z • {^^(pMp) + P^^HpMp)) ■ ■ (2.70) 

p 


A difficulty with this interpretation is that the operator ^(0) is not defined. This does 
not matter for Kronig’s identity, but it is a problem for the CCR. One natural solution 
to this is as follows. The non-trivial commutation relation in (2.35) implies 

(R+R_)(Q+ + Q_)(R+R_)t = g+ + g_, (2.71a) 

{R+R-)U{0){R+R-)^ = fl(0) + V2. (2.71b) 

This suggests to formally identify R+R- with 

( 2^/2ti 

exp -i—f— 


1 .( 0 ) 


(2.72) 


since this is compatible with (2.71) and the CCR in (2.69) for p = p' = {). The natural 
interpretation of d*(0) is thus in the form of a variable on a circle with radius L/{2^/2tt), 
i.e. $(0) should be identified with ^>(0) + Lj\/2. We note that this interpretation is also 
compatible with the discrete spectrum of 11(0). 


3 From bosons to fermions 

This section is dedicated to the reconstruction of the fermion held operators (k) from 
the fermion densities Jr{p) and the Klein factors Rr- The main objects of study will be 
so-called vertex operators 

xR)" exp —a{p)Jr{—p^ x (3-1) 

^This Hamiltonian corresponds to an important model in particle physics describing relativistic bosons 
in 1-1-1 dimensions; cf. Remark B.l in Appendix B. 
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with integers w and complex parameters a{p). The crosses in (3.1) indicate boson normal 
ordering which is defined in (3.6) below. Our goal is to prove that the vertex operators 

Pl^-ipx-e\p\/2^ j X ^3 2) 

converge in the limit e —)■ 0"*“ to the formal fields iprix) in the introduction, in a sense 
made precise in Proposition 3.7. Using this we will also develop tools needed to com¬ 
pute fermion correlation functions. We note that the formula in (3.2) appeared in the 
condensed matter physics literature in [23, 26] and in the mathematics literature in [34]. 

Remark 3.1. The name of the objects in (3.1) is motivated by the study of vertex operator 
algebras in the mathematics literature; see e.g. [20]. 0 


A{x]e) = 


Vl 


kR^ exp 


27rf 


jr{o)x — 

p¥=o 


Jr{- 

ip 


3.1 Vertex operators 

We consider periodic functions on the space [—L/2, L/2] given by 

m = w‘^y + ^a(p)e*™ {y G [-L/2, L/2]), (3.3) 

p 

where w € Z and a{p) G C satisfying the following conditions; 

_ _ <2^ 

a{p) = a{-p), y]] —|p||a(p)|^ < oo. (3.4) 

p 

The first condition is equivalent to requiring that / must be real-valued, i.e. f{y) = f{y) 
for all y G [—L/2, L/2], and the second is a weak smoothness condition on /.^ We denote 
by Q the set of all functions satisfying these conditions. 

For each G we define 


= ^iTra(0)JriO)/Lj^w^ina(0)Jr(0)/L 


\ PT^O 



(3.5) 


If / is real-valued, then the properties of the fermion densities Jr{p) and the Klein factors 
Rr derived in Section 2 suggest that r(e*'^) defines a unitary operator. As we will show, 
this is indeed the case. 

Remark 3.2. We note that, in order for e*-^ to be periodic with period L, it is enough 
to require that / is periodic on [—L/2, L/2] up to an integer multiple of 27r. It follows 
that, since w = [/(L/2) — /(—L/2)]/(27r) and since a{p) are the Fourier coefficients of 
the periodic function f{y) — w27rx/L, the parameters {w, {a{p))p) are in a one-to-one 
correspondence with such functions e*-^. This integer w is often referred to as winding 
number. 0 

^Functions satisfying the second condition in (3.4) are sometimes said to be this condition is 

weaker than differentiability, but stronger than continuity. 
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Remark 3.3. It is easy to verify that Q forms a group under multiplication (with ^Re^’R = 
g*[/i+/ 2 ]^ etc.). In fact, this group is an important example of a so-called 

loop group, and Tr{e^-^) in (3.5) corresponds to a projective representation of Q [5]. 0 

The vertex operators in (3.1) are related to Tr{e^-^) in the sense that the former can 
be written as xrr(e*-^)x . The precise definition of x - • - x for such operators is 

xr^(e*^)x = (3_g) 


with a short for {a{p))p and 

J^{a) = ^ ‘^a{±rp)Jr{^rp), j)?(a) = ^a(0)J^(0). (3.7) 

p>0 

As a motivation, we note that Proposition 2.5 suggests that^ 

J~ (a)n = (a)^II = 0, (3-8) 


and it is therefore natural to regard (a) and J~ {a) as creation and annihilation 
operators, respectively. We conclude that (3.6) indeed is a normal ordering prescription; 
all creation operators are to the left of the annihilation operators, and the zero-mode 
parts, namely J)?(q;) and the Klein factors Rr, are placed symmetrically in the middle. 

We find it convenient to introduce the shorthand notation 


Cr{ai,a2) 


^ ‘^pai{-rp)a2{rp), 

p>0 


g*/2) = gi7rr(«;2«l(0)-t«ia2(0))/Lg-Cr(ai,O2) 


(3.9) 

(3.10) 


Note that the second condition in (3.4) is equivalent to requiring that Cr{a,a) is finite. 
Moreover, if both Q;i(p) and 02 (p) satisfy this condition, then Cr{ai, 02 ) is finite (by the 
Cauchy-Schwarz inequality). 

The following proposition summarizes important properties of the vertex operators 
which we will need: 

Proposition 3.4. For each G Q, Tr{e''f) in (3.5) defines a unitary operator on T 
satisfying 


r,(e*^) = ^ 

(3.11a) 

X X 

T 

[-H 

X X 

II 

X X 

Ph 

X X 

(3.11b) 


(3.12a) 

i?±xr^(e*-^)x = (-l)"'xr^(e*-^)xi?±, 

(3.12b) 


®As will be shown, the jf{a) are well-defined as quadratic forms on T>p, and this gives a precise 
meaning to the relations in (3.8); see Section 3.4 and also Appendix C. 
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together with the following commutation relation on Vp: 

Jr'{p),iTr{e^^)i = rSr^ {w5pfi + ipa{p))xTr{e''^)x . (3.13) 

Moreover, the following identities hold true: 

xr±(e*^)xxr±(e*^')x = C±(e*-^,e*-^')xr±(e*[-^+-^'l)x , (3.14a) 

xr±(e*^)xxr^(e*^')x = (-l)^’^'xr^(e*^')xxr±(e*-^)x , (3.14b) 

and 

(ll,xr+(e*-^)xxr_(e*-^')xll) = 5^,o5«;',o (3.14c) 

for all e*-^, e*-^' G Q. 

(The proof is given in Section 3.4.) 

Eqnation (3.11a) implies that the operator xrr(e*-^)x is bounded for all e*'^ G Q (since 
it is equal to a unitary operator multiplied by a finite constant). Therefore, as long as 
the conditions in (3.4) are satisfied, formal computations with vertex operators are safe. 
However, as we will see, in order to reconstruct the fermion field operators using vertex 
operators, it is necessary to pass to the non-trivial limit \a{p) \ —)■ l/|p|. 

We now state a simple implication of (3.14a), which, together with (3.14b)-(3.14c), 
functions as a working horse for computing fermion correlation functions. 

Corollary 3.5. Let iV > 1 6e on integer and consider e^^i G Q for j = 1,2,... , N . Then 
the following holds true: 


<r^(e*^i)x...xr^(e*^^)x = 


with Cr{-,-) defined in (3.10). 


Yl a(e*^^e*^i') xr,(e*[^i+-+^^l)x (3.15) 


(The proof follows by a simple induction argument.) 

As a closing remark we mention that, for the formal fermion field operators 'iprix) in 
the introduction, the following identity holds true: 


rr{e^^)^‘f’r'ix)Tr{e^^) = 


if r' = r, 
'ipr' ( 2 :) if r' 7 ^ r 


(3.16) 


(this follows from (3.2) using identities given in Proposition 3.4). The r.h.s. of this has 
a natural interpretation as a chiral gauge transformation, and, for this reason, Tr{e^^) 
are sometimes referred to as implementers of gauge transformations in the literature. 
In the case of the space [—L/2,L/2], there exist so-called large gauge transformations 
where /(L/2) — /(—L/2) is a non-zero integer multiple of 27r. It follows that the Klein 
factors Rf correspond to the implementers of large gauge transformations given by the 
functions f{y) = w2TTy/L. It is also interesting to note that (3.13) implies 

Tr{e^^)^ Jr'{x)Tr{e^^) = Jr'{x) + rd^y —dxf{x) 


(3.17) 
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for the operator-valued distributions Jr{x) related to jr{p) by the Fourier transform in 
(B.12) (cf. also (B.3)). 

Remark 3.6. The observation that the Klein factors Rr and the fermion densities Jr{p) 
function as implementers of gauge transformations makes clear that they should be 
treated on an equal footing. This is also the reason why we sometimes use the word 
bosons to refer to both of these operators together. 0 


3.2 Reconstruction of the fermion field operators 


Let e > 0 be a regularization parameter and consider functions 

Ueiy) = - ^) + E (y G [-L/2, L/2]) (3.18) 

^ P^o ^ 

parametrized by a variable x G [—L/2, L/2]. Since these functions satisfy the conditions 
in (3.4), it follows that they can be used to define vertex operators . Our 

main result in this section is that these allow us to reconstruct the fermion field operators. 

Proposition 3.7. The following holds true as an identity on Vp-f' 

= dx (3.19) 

(The proof is given in Section 3.5.) 

One important observation is that the r.h.s. in (3.19) is a well-defined operator map¬ 
ping T)p to Vp (since this is true for the l.h.s.). Combining this with the Fourier trans¬ 
form in (B.12) justifies (3.2) since the r.h.s. in the latter is equal to (1/\/L) ^ 

and it is natural to denote this operator by f^rix'-, e) since it converges to f^rix) in the 
limit e —7- O"*" in a sense made precise by (3.19) and Remark 3.9. 

For the sake of computations, it is useful to write (3.2) in the following form: 


V'r(x;e) = 

V27r€ 


r E (3.20) 

I n 




Lp 


with 




1 _ g-Lc/i - 1 + 0{e/L) 


(3.21) 


(this follows from (3.11a); cf. also (A.l)). Since one is usually only interested in the limit 
e —7- 0, the constant A// is often ignored. For later reference we also state the Hermitian 
conjugate of (3.19): 

fL/2 

(3.22) 

I-L/2 


1 

iblik) = lim , / dx 

^^oV^J-L/2 

(this follows from (3.11b)). 


The precise meaning of the limit in (3.19) is given in Remark 3.9. 
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Remark 3.8. Note that fx,e{y) are C'°°-functions which converge to the distribution 
7rsgn(y — x) in the limit e —)• 0^. This is the key observation which allow the operator¬ 
valued distribution ipr{x) to be obtained as a limit of a bounded operator iprix] e), which, 
as we will see later, is very useful to define and compute fermion correlation functions. 

0 

Remark 3.9. The limit in (3.19) is in the strong operator topology [30]. We emphasize 
that one can use other regularizations than the damping factor exp(—e|p|/2) in (3.18). 
For instance, this factor can be replaced with a sharp cutoff 0(A — jpj). This would have 
the advantage that the corresponding limit A —)■ oo in (3.19) would be in the sense of 
Definition C.5, which is somewhat stronger. However, the e-regularization is easier from 
a computational point of view. 0 


3.3 Correlation functions 


To introduce some notation and techniques needed in Section 4, we demonstrate how 
the results in Sections 3.1 and 3.2 can be used to compute fermion correlation functions. 
For simplicity, we restrict ourselves to equal-time correlation functions. 

Proposition 3.7 and the discussion directly thereafter suggest to use the regularized 
fermion operators ^jJr\x;e) in (3.2) to define the correlation functions 

(D,V’^((xi;0+).. = lim (D, V’?) (xi; ei)... V’?)^(xAr; eAr)H) (3.23) 

for rjjQj = ± and xj G [—L/2,L/2], where we use the following convenient notation: 

ilj+{x-,e) = i^l{x;e), i^~{x;e) = i;^{x;e). (3.24) 


As will be seen, with this definition, one automatically obtains prescriptions of how to 
interpret certain singular functions as distributions; see Proposition 3.10. Moreover, it 
follows from (3.19) and (3.22) that V’r(x;e) = with fx,e{y) in (3.18). 

To compute the fermion correlation functions we use the following two special cases 
of (3.14a) and (3.14c), respectively: 


V’^(x;e)V^«'(x';e') = 


2 sin 


iexp(^(e + eO) 
[l[r{x-x') + \{e + e% 


-qq 


x'il)‘^{x]e)'il)‘j. {x']e')x (3.25) 


(this follows using (A.l)) and 

(H, (xi; 0+)... 0+)(Al) = Rl\^^ ... (3.26) 

The r.h.s. in (3.26) is a sign factor which is either 0, 1, or —1, and we note that it is zero 
unless both qnSr„,+ and J2n=i QnSr„- are zero (this follows from the identities in 

Corollary 2.8). In particular, 




^ri,r2^qi,—q2 ) 

^ri,r2^r3,r4Sqi,-q2^q3,-q4, ~ ^ri,r3^r2,r4^<}l,-<l3^g2,-q4 


(3.27) 

(3.28) 


Using Corollary 3.5 we get the following result: 
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Proposition 3.10. The fermion correlation functions in (3.23)-(3.24) are equal to 

/ , \ ~Qnl3m^rn,rm 


(Q, TT - , , ' -r-— 

l<niLiV J 

In the proof of this, and also in Section 4, we need the following: 


(3.29) 


?giri 


zero for some N gN. Then 


Q.nQin^rn,rm ~ ' 

l<n<m<N l<n<m<N 


thus E„<m QnqmSr„,rm is equal to 

r n \rn^n j 


Qnhr„,r ( Qnhr„,r) 


Similarly, En<m = (l/2)En,r9n5r„,r (-g'n'5r„,-r) = 0. 

Proof of Proposition 3.10. We are left to show that 



is non- 

= 0. 

(3.30) 

0 for r = 

±, and 

N 

(3.31) 

T‘ 


□ 


N 

INL\ = 


iqnqmSr„,r^ if ^ ^ 0, 


l<n<m<N 

but this is true according to Lemma 3.11. 


(3.32) 


□ 


Note that, due to (3.32), the thermodynamic limit L —)■ oo of the correlation functions 
in (3.29) exists and is equal to 




1 


n 


l<n<m<N 


ri^Xfi Xyo) “t“ i0~^ 


qnqmSr„,rm 


(3.33) 


Remark 3.12. A careful reader might wonder how the result in Proposition 3.10 can be 
consistent with Wick’s theorem, which applies to the present model and states that any 
A^-point correlation function can be expressed as a linear combination of products of 
two-point correlation functions. It can be shown that one only needs to check the cases 
N = 2M with (rn, Qn) = (r, -|-) and {rM+n, QM+n) = (r, —) for 1 < n < M . The result in 
Proposition 3.10 is then equivalent to the one obtained by Wick’s theorem if and only if 


Ul<n<m<M sin(t/n “ Um) sin(F^ - Vn) 


Un,m=l Sin(t/n - Vm) 


= det 


1 


l<n,m<M \sin(Ln ” Im) 


(3.34) 


for Un = vrXn/L and Vn = nx 2 M+i-n/L. The latter holds true due to a well-known 
variant of Cauchy’s determinant identity. 0 
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3.4 Proof of Proposition 3.4 


The identities in Proposition 3.4 can easily be derived using the following special cases 
of the Baker-Campbell-Hausdorff formula:^ 


R 


if [A, B] = cl, c G C, 


(3.35) 


[A, e^] = ce^ 


together with properties of the fermion densities and the Klein factors summarized in 
Proposition 2.5 and Corollary 2.8; see e.g. [41]. However, this would not constitute a 
proof; the J^{a) are unbounded operators, but the identities in (3.35) can only be used 
in such a way for bounded operators. Moreover, it is not obvious how the exponentials 
of the unbounded operators J^{a) should be defined (or that they can be defined at all). 
To prove these identities we give a precise meaning to the vertex operators in (3.6) as 
quadratic forms; see Appendix C. We first concentrate on vertex operators with trivial 
zero-mode parts, i.e. q;( 0) = tc = 0. 

Our strategy of proof is as follows. We start by showing that the vertex operators 


XgiJ+(o)+iJr (o)x _ (o) 

are well-defined as quadratic forms onDp and that the identities 

(x^iJ+(o)+ijy(a)x)t _ x^-iJ+(a)-iJ-{a)^ ^ 

X iJ+(ai)-|-iJ“(ai)x X ij+ (a2)+iJr ( 012 )^ _ -Cr(oi,02) ^piJr iai+a2)+iJr (ai+a2)> 

xo xxo xo xc y 


(3.36) 

(3.37a) 

(3.37b) 

(3.37c) 


hold true for such quadratic forms; see Lemma 3.14. Note that, formally, the definition 
of J^{a) in (3.7) and the commutation relation in (2.21c) imply 


Jr {oil), {oi2)\^ = Cr{oil,a2), 

Jr{p), Jr (a) + Jr{^) \ = irpa{p), 


(3.38a) 

(3.38b) 


with Cr{ai,a 2 ) in (3.9). The identities in (3.37b) and (3.37c) thus follow from the first 
and third identities in (3.35). However, we need to justify the use of these identities for 
quadratic forms. 

In the next step, we show that, for a = {a{p))p satisfying the conditions in (3.4), the 
quadratic form 

{a)+ijp {a) — ^-Cr{a,a)/2 {a)+ijp (a)^ (3.39) 

can be identified with a unitary operator on T"; see Lemma 3.17. (We emphasize that we 
regard (3.39) as a definition motivated by the second identity in (3.35).) Since there is no 
risk of confusion, we can thereafter use the same symbol (a) fQj. unitary 

operator. With these results in place, the proof of Proposition 3.4 is straightforward. 


^See e.g. Appendix C in [41] for proofs of these identities valid for bounded operators A and B. 
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Remark 3.13. It is interesting to note that our proof of (3.37a)-(3.37c) does not require 
the conditions in (3.4): the quadratic form in (3.36) is well-defined for any a = {a{p))p, 
i.e. no reality or decay condition is needed, and (3.37a) (with a on the r.h.s. replaced by 
a) and (3.37c) are always true. Moreover, (3.37b) holds true provided that 

'^p\ai{-rp)\\a 2 {rp)\ < oo. (3.40) 

p>0 

While we do not make use of these more general results, they are important in conformal 
field theory; see e.g. [20]. 0 

A simple way to make the vertex operators in (3.36) precise is as quadratic forms, 
with the exponentials defined as Taylor series: 

Lemma 3 . 14 . Let a = {a{p))p satisfy the conditions in (3.4). Then the vertex operator 

X iJ+(a)-hiJy(Q)x ^ 'Sp [iJr{(^)T' [iJL («)]"" ('qzLIi 

^ - 2^ nl m\ ^ ’ 

n,m=0 

is well-defined as a quadratic form on Vp. Moreover, for ai and 02 similarly defined, 
(3.37a)-(3.37c) hold true as identities of quadratic forms on Vp. 

Proof. Note that, for boson states and as in (2.56), 

(,®.(j+(a)r(jr(a)r*> = {(jr( 5 ))”i?S.(.'r(a)r»S-) (3.42) 

is well-defined for all (n, m) G Nq and that it is non-zero only for finitely many such 
pairs (this follows from (2.21c) and (3.8)). This together with Proposition 2.12 implies 
that (3.41) is well-defined as a quadratic form on Vp. 

As a consequence of this and = Jfi{a) (this follows from (2.19) and (3.7)) 

we conclude that (3.37a) holds true. 

To prove (3.37b), we note that, by a similar argument as above, 

yrjGVp. (3.43) 

Indeed, for all boson states as in (2.56) and n G N, {Jfi is a linear combination 

of such boson states, and it is non-zero only for finitely many n (this follows again from 
(2.21c) and (3.8)). Thus, by Proposition 2.12, /n\){iJfi (a))"'?] is well- 

defined and in Vp. Therefore, to prove the validity of (3.37b) in the sense of quadratic 
forms, we only need to show that 

{ri,e^e^ri') = e^{'q,e^e^fi) y-i],?]' eVp (3.44) 

for A = iJfi (cci), B = iJfi{a 2 ), and c = —Cr{ai, 02 )- We note that (3.38a) is equivalent 
to [A, S] = c, which is obviously well-defined as an identity of quadratic forms on Vp. 
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Moreover, this implies that® 


min(n,m) 

='^11 

1=0 

holds true as an identity of quadratic forms on Bp. Recalling that (ry, B^A^r]') is non¬ 
zero only for finitely many (n,m) G N^, one can expand the exponentials on both sides 
in (3.44) as Taylor series and use (3.45) to compute the l.h.s. One thus finds that both 
sides are well-defined and identical.® 

The proof of (3.37c) is similar, but simpler, and thus omitted. □ 


(3.45) 


Lemma 3.15. Given a = {a{p))p satisfying the conditions in (3.4), the quadratic form 
^iJf{a)+iJr (a) (3.39) can be identified with a unitary operator on T. 


Proof. Using the definition and (3.37a)-(3.37b) we compute 


ijf (a) _ ^—Cr{a,a) J^(Q;)-|-i J,. (a)j 


x6 


' xxC 


= L 


and similarly 


(3.46) 


^ijf {a)+iJr (o)^ ^ijf{a)+iJr («) _ J |'3 

This and a general result for quadratic forms (see Lemma C.4) conclude the proof. □ 


Proof of Proposition 3.f. It follows from the self-adjointness of Qr and the unitarity of 
Rr that (a)/2 unitary operator. Equations (2.21c) and (2.38) imply 

that this operator commutes with («)^ From this and Lemma 3.15 it follows 

that F^(e*'^) in (3.5) defines a unitary operator on 7^, and (3.39) then implies that (3.11a) 
holds true. The identity in (3.11b) thus follows directly from (3.4)-(3.5) and (3.37a). 

The remaining identities can be proven using Proposition 2.5, Proposition 2.7, and 
Corollary 2.8. As an example we prove the last three identities. Equation (2.41c) implies 
that 


giJ°(ai)/2^«)lgij“(oi)/2gij0(a2)/2^ui2giJ°(«2)/2 

_ gi7rr-(«)2Ol(0)--!«ia2(0))/LgiJ°(ai+O2)/2^wi+-!«2giJr (oi+a2)/2 


®The identity in (3.45) can be proved by induction over m £ No for fixed n £ No using standard 
relations for the binomial coefficients. 

®The key identity is 


EE= E 


N M min(Al,M) 


^ c{y,B"'A"y') {N,M€N). 


n =0 m =0 1=0 


n\m\l\ 
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which together with (3.37b) yields the identity in (3.14a). The identity in (3.14b) follows 
from (2.41a). Lastly, we note that the l.h.s. of (3.14c) is equal to 

= (fl,e*‘^+(")/2gijO(a')/2^»^»'giJ^(a)/2giJ“(a')/2^^ = {VL,R'IR^VL) (3.49) 
(where we have used (3.8) and Qr^ = 0), and the r.h.s. then follows from (2.42c). □ 

3.5 Proof of Proposition 3.7 

We divide the proof into three steps. In the first step we show that the r.h.s. in (3.19) 
is a well-defined operator, which we denote by Vr{k), satisfying a list of four identities 
(see Lemma 3.16) which are identical with identities satisfied by the fermion field oper¬ 
ators 'iprik)- In the second step we show that all matrix elements of Vr{k) with respect 
to fermion states in (2.6) can be computed using only the identities in this list (see 
Lemma 3.17). We thus conclude in the third step that both sides in (3.19) agree as 
quadratic forms on Vp, and we argue that this is only possible if they agree as operators 
on Dp. 


Lemma 3.16. Denote the r.h.s. in (3.19) by Vr(k). Then Vr{k) is a well-defined operator 
mapping Vp to Vp and satisfying the following identities: 


Vfirk)D = V}{-rk)D = 0 VA: > 0, 

(3.50) 

with Vr {k) = Vr{k)\ together with 


Jr ip) > (^) = ^r,r' V^{k-p), 

(3.51) 

R±V±{k) = V±{k + ^)R±, R±Vp{k) = -Vp{k)R± 

(3.52) 

on Vp, and 


{Rln,v4-l)n) = yj, {R_n,v_{i)n) = 

(3.53) 

(The proof is given below.) 

It is important to note that the identities in Lemma 3.16 are also satisfied by the 
fermion field operators fir{k). (3.50), (3.51), and (3.52) correspond to (2.3), (2.21a), and 
(2.43), respectively, and the identities corresponding to (3.53) are 


(3.54) 


(the latter follow from (2.2)-(2.3) and (2.43)-(2.44)). 
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Proof of Lemma 3.16. The vertex operators xr^(e can be written in the form 


F±(x; e) = n exp J,(^rp)eW-ep/ 2 ^ . 


(3.55) 

(3.56) 


We note that F^’^(x;e) and Vp{x;e) correspond to creation and annihilation operators, 
respectively, since 

V~ {x-,e)Ll = V^{x-.,e)^Ll = Ll, (3.57) 

as suggested by (3.8) (this is proved below). The exponentials in (3.56) are well-defined 
as Taylor series: we can represent ^^^(x; e) as sums over integer vectors u± = {i'±{p))p>o 
with elements v±{p) G Nq, 


V,^{x-,e) = Y,Ur{u±)e^^' 


irxP(u±)—eP(u±)/2 


(3.58) 


using the shorthand notation 




p{^±) = '^p^±{p)- (3.59) 

p>0 


We note that (3.58) is well-defined as a quadratic form on for fixed finite linear 
combinations p and p' of boson states in (2.56), (p,xrr(e~^^P’^)x p') can be computed 
using (3.55) and (3.58)-(3.59), and one finds that only a finite number of the terms 
in the resulting i/±-sums are non-zero, which implies that the same holds true for all 
p,p' G Vp (cf. Proposition 2.12). Moreover, using (3.58)-(3.59), the identity in (3.57) 
follows from Proposition 2.5. 

We now fix a boson state 


pg = B^RfRj-n, = n 


5t(p)m(p) 


(3.60) 


and compute xrr(e ^'^P’^)xp^ with the help of (3.55)-(3.59). It follows from (2.41c) 
that 

p^Q-RfRZ'^-n = e^‘^‘^^RfRZ‘^-n VaGM, (3.61) 

which together with the above identities yields 




with the shorthand notation 


Wriqr, {i^±}) =‘^(qr-^'^ -r [P{up) - P(i/-)] . 


(3.63) 
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The integral over [—L/2,L/2] in (3.19) and the i^^-sums can now be exchanged (this is 
allowed due to the e-regularization). Hence, for fixed k, the integral 

1 / r 

_ ^ (3.64) 

imposes the following constraint on the iz±-sums in (3.62): 

k=‘^ (qr-^) -r[P{u+) - P{u.)]. (3.65) 

For any boson state, there is at most a finite number of iz_-terms which give a non-zero 
contribution to (3.62) (this follows from (2.20) and (2.56)). Hence, there exists some 
P < oo such that P(iz_) < P for all iz_-terms. It follows that (3.65) restricts the iz_(_- 
sum in (3.62) to a finite number, which means that all non-zero contributions are given 
by pairs (i/_|_,iz_) belonging to some finite set determined by the constraint in (3.65). 

In conclusion, both iz±-sums in (3.62) are finite, and the limit e —)■ 0^ in (3.19) can 
therefore be taken without problem. The result is 

yr{k)v^ = (3-66) 

,Z/_ 

where the tilde above the summation symbol indicates that the iz^-sums are restricted 
by (3.65); the corresponding formula for Vr {k)ri^ can be derived in a similar manner 
(we omit the details). It follows that Vr{k) is a well-defined operator mapping Vp to 
Pp, and we are only left to show that this operator satisfies the stated identities. 

Let r/® = n in (3.60) and use (3.66) to compute Vr{k)Q,. All non-zero i/_-vectors will 
give a zero contribution to the iz_-sum in (3.66), and the constraint in (3.65) therefore 
simplifies to 

k =—r^ — rP{up). (3.67) 

1j 

If rk > 0, then no iz+-vector can satisfy (3.67), which implies that 

Vr{h)P = 0 yk : rk > 0. (3.68) 

This proves the first identity in (3.50); the second follows from a similar argument applied 
to Vr {k)^. li k = —rp/L, only the solution P{up) = 0 is allowed, in which case (3.66) 
yields 

Pip - Pip - 

(3.69) 

This and the unitarity of imply (3.53). 

Equation (3.13) implies that 

Jr{p)PTrPe-^^'R-p^ = , (3.70) 
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which together with (3.19) yields 


Jr (p)) ^r' iJJ) 


—6ryVr>{k + p). 


(3.71) 


By taking the adjoint of this identity and using (2.19) we obtain (3.51). 

From (3.12a) and (3.12b) we obtain 

77±xr±(e^*-^"'')x = (3.72a) 

i?±xr^(e^*-^"’')x =(3.72b) 

and this yields (3.52) by recalling that Vr{k) is defined as the r.h.s. of (3.19). □ 


Lemma 3.17. For all rj, r/' G Vp, the matrix element [r], Vr{k)rj') is uniquely determined 
by the identities in (3.50)-(3.53). 

Proof. Due to Proposition 2.12 we can restrict ourselves to the boson states in (2.56). 
We thus set rj = and rj' = which we write as (cf. (3.60)) 


Vm = ^(m,+)^(m-)7?+ Rl D, 


p>0 


Ijt (^J.p^m{rp) 

\/m{rp)\ 


(3.73) 


and similarly for rj^,. The matrix elements {r]^, Vr{k)r]^,) can thus be written 


which is non-zero only if q_^ = q'_^ and qr = <l'r ~ (this follows from (2.41a), (2.42c), 
and (3.66) since R± and 6^(p) commute). Furthermore, from (2.32) and (3.71) it follows 


that 


Vr{k)b\r'p) = b\r' p)Vr{k) + irSry 



rp), 


(3.75) 


which implies that the matrix elements in (3.74) are non-zero only if i?(ni,-r) = R(m',-r)- 
It follows that we can restrict ourselves to states rj^ and rj^, with {q^, q_^) = {q—r, 0) 
and {q'r-,q'-r) = (9)0) for some g G Z and 71(m,-r) = R{m',-r) = 7: all non-zero matrix 
elements in (3.74) can be written in the form 


{B^m,r)Rr"^, Mk)B(m',r)Rt^)- (3.76) 

By repeated use of (3.75) we can write (3.76) as a linear combination of matrix elements 

{k - rP(m")) (3.77) 

where P(m") is defined as in (3.59) and m" is any integer vector satisfying 0 < mf{p) < 
m'{rp) for p > 0 and mf{p) = 0 for p < 0. Using (2.33), we see that (3.77) is non-zero 
only if = B(^^i_^n j.y Hence, only the matrix elements 

Vr {k - rP(m")) R^r^) 


(3.78) 
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need to be considered. These can be computed using (3.52) as follows: 


Vr (k - rP(m';)) R^n) = (R-^Q, R-^Vr {k - rP(m")) R^r^) 

= {R;''n,Vrik')n), 


(3.79) 


where k' = k — — rP(m('). We are left to compute the matrix element on the r.h.s. 

of (3.79). For rk' > 0, this matrix element is zero due to (3.68). For rk' = —TxjL^ it 
is equal to ^jLj (27r) due to (3.53) and the unitarity of Rr- For rk! < —tt/L, we can 
rearrange the r.h.s. of (3.79) as 

(i?-’'f2, Vrik!)^) = (O, R:Vr{k')R-^Rl^) 

= {^,Vr{k' + r^)Rl^) = {V}{k' + r^)VL,Rin) 

(the second identity follows from (3.52)), and it thus follows from (3.50) that this is zero 
since r{k' + r^) <0. □ 

Proof of Proposition 3.7. Lemmas 3.16 and 3.17 imply that 

(? 7 , Vr{k)ri') = {rj, 'tpr{k)r]') Vr?, r?' G Vp, (3.81) 


and by general results for quadratic forms it follows that (3.19) holds true as an operator 
identity on Vp (see Appendix C). This concludes the proof of Proposition 3.7. □ 


4 Fermion-phonon model 

In this section we use the mathematical results developed in Sections 2 and 3 to study 
the fermion-phonon model described in the introduction. As will be shown, this model 
is exactly solvable in the sense that not only the eigenstates and eigenvalues of the 
Hamiltonian, but also all correlation functions, can be computed by analytical means; 
we will, however, be content with discussing fermion correlation functions. 

4.1 Definition and snmmary of resnlts 

We give a precise definition of the fermion-phonon model in Fourier space using fermion 
field operators 'ipi^\k) (r = ±) and phonon field operators 4>p(p) and IIp(p). We recall 
that L > 0 is an infrared cutoff which constrain the fermion and boson momenta such 
that k G {2pjL){'Z + 1/2) and p G {2pfL)'Z, respectively. The interactions are given by 
the fermion-fermion and the fermion-phonon coupling constants A and g, respectively, 
and they are regularized using an ultraviolet cutoff a > 0; see (4.3) below. The other two 
model parameters vp > 0 and up > 0 with vp < vp correspond to the Fermi velocity 
and the phonon velocity, respectively. As will be shown, the four model parameters must 
satisfy the conditions in (1.4). For reasons explained in Remark 4.3, we also introduce a 
less important infrared cutoff ujq > 0. We remind the reader that some of our notation 
is summarized in Appendix A. 
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Definition 4.1. The fermion-phonon model is defined by the Hamiltonian 


H = vfJ 2Y1 ■i’l{k)'4>r{k): + X] 7 ' + <^p{pf^p{p)^ 


r k 


L 


E X (k^JliPUp )+E ^^>pJHp)ip{p)] (4.1) 


^ ^ L \ 27r 

p \ 


with the fermion densities Jr{p) = ■'tpr{k)il’r{k + p):, where the colons indicate 

normal ordering defined in (2.10), the free-phonon dispersion relations 


u%{p) = 


cjo if p = 0, 

vp\p\ if P / 0 , 


and the regularized interaction potentials 

Xa{p) = A0(7r/a - IpI), paip) = gQ{Tr/a - \p\). 

The fermion and phonon field operators satisfy the following CAR and CCR: 

r 7 ... 7 + .,,.'1 T . . 




L 


where V'J(/c) = '4>r{k)\ 


[^,.(fc),^p(/c')^ 

•H 

^l{k)Al{k')\ 

^p{p)AP{p) 

= 

np(p),rip(p')' 

^p{p),^i^\k) 

= 

np(p),7(^)(fe)' 


(4.2) 


(4.3) 


= 0, (4.4c) 


<f>f,(p) = l>p(p)f = 4>p(-p), n^(p) = UpipY = Up{-p) 


t = 


_ 


p\-p)^ 


(4.5) 


and these are represented on a Hilbert space J- which contains a normalized state 14, i.e. 
(n, 14) = 1 with the inner product (•, •) in 7, such that 

'ip^{rk)Vt = ^l{—rk)Vt = 0 V/c > 0, 

/ . „ . X (4.6) 

(^inp(p) -h Wp(p)4>p(p) j 14 = 0 Vp. 

We refer to F as the fermion-phonon Fock space. 

It is important to note that we use the regularized interactions in (4.3) only as a 
specific example to simplify our notation. As will become clear, all computations in the 
regularized model can be done without specifying the explicit form of Xa{p) and ga{p)- 
In fact, our results for the regularized model hold true for a large class of such potentials; 
the only restriction is that these potentials vanish sufficiently fast as |p| —)■ oo and that 
|Aa(p)| < |A| for all p and similarly for ga{p)- Moreover, we expect that our results are 
universal in the sense that they are (essentially) independent of the regularization in the 
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limit a —)■ 0+ as long as lim^_^o+ ^a{p) = A (independent of p) and similarly for ga{p)- 
However, a proof of this is beyond the scope of the present paper. 

We begin with a summary of the results obtained in this section. The first set of 
results concerns the regularized model with cjq > 0, a > 0, and L < oo: 

• We show that the fermion-phonon Hilbert space J- is fully specified by the condi¬ 
tions in (4.4)-(4.6); see Section 4.2. 

• We prove that the Hamiltonian H in (4.1)-(4.3) defines a self-adjoint operator on 
T" with pure point spectrum and non-degenerate ground state, and we construct 
all its eigenstates and eigenvalues; see Corollary 4.6. 

• We develop tools by which any fermion correlation function can be computed; see 
(4.33)-(4.36). 

The second set of results concerns the limit where cuq —?• 0"*“, a —)• 0"*“, and L —)■ oo, which 
leads to simpler formulas. The physical motivation is that this limiting description is 
adequate for describing properties at intermediate length scales which are much larger 
than the interaction range a and much smaller than the system size L. The price to 
pay for this simplification are divergencies which appear in formal treatments. However, 
these divergencies can be eliminated by additive and multiplicative renormalizations: 

• We construct a renormalized fermion-phonon model in the continuum limit a —)• 0"*“; 
see Proposition 4.8. Our construction makes manifest that the Hilbert spaces for 
a > 0 and a = 0 are not unitarily equivalent. 

• We obtain explicit formulas for all fermion correlation functions in the limit a —)■ O"*" 
and L —)■ oo; see Corollary 4.10. 

4.2 Fermion-phonon Fock space 

To construct the Hilbert space J- from the vacuum H we follow the method in Section 2. 
In the present case, we need both fermion and boson creation operators ct{k) and 6p(p), 
respectively; these are operators on T satisfying the conditions 


(fc)) C^/(fc ) ^ — br,r'bk,k'j 

|c^(A:),Cp(/c')| = 0, Cr{k)n = 0, 


bp{p),b^p{p) =5py, 

bpip),bp{p) = 0 , bp{p)n = 0 , 

(4.7) 

Crik),bp\p) =0, 

clik),b^p\p) =0, 



with Cr{k) = cl{k)^ and bp{p) = 6p(p)^. Using these we can construct the following 
states in T: 

(Jl G {0,1}, mp(p)G No, (4.8) 
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where only finitely many of the quantum numbers nr^F{k) and mp{p) are non-zero, and 
where n is short for the infinite vector {{nr^F{k))r,k, irnp{p))p)-^^ The relations above 
imply that the states in (4.8) are orthonormal, and they also determine the action of the 
creation and annihilation operators on these states. Let T> denote the vector space of all 
finite linear combinations of the states in (4.8). This set forms a pre-Hilbert space, and 
T" is obtained from D by norm-completion [30]. 

Remark 4.2. The Hilbert space J- thus constructed can (of course) be identified with the 
tensor product Tf <8> Tp of a fermion Fock space and a phonon Fock space. The Fock 
space Tf is identical with the fermion Fock space constructed in Section 2.1, and the 
states in (4.8) correspond to a tensor product of the fermion states in (2.6) with phonon 
states. It is thus clear that all bosonization results developed in Sections 2 and 3 can be 
used for the fermion-phonon model. 0 


The representation of the field operators used in our definition of the fermion-phonon 
model is such that the free fermion-phonon Hamiltonian iLg, which is obtained from 
the Hamiltonian in (4.1) by setting Aa(p) = Qaip) = 0, is a (positive-definite) self- 
adjoint operator with eigenstates in (4.8). To see this we note that the fermion field 
operators can be constructed as in (2.9), and for the phonon operators the corresponding 
formulas are 


^p{p) = -i 




^p{p) 



(4.9) 


Indeed, one can verify that this implies that all the conditions in (4.4)-(4.5) are fulfilled 
and that the free fermion-phonon Hamiltonian can be written as 


Ho = ^uf|/c|cJ(/c)c^(/c) + J2‘^piP)^p(p)bp{p)- 

r,k p 


(4.10) 


Moreover, as the latter formula makes manifest, the states in (4.8) are exact eigenstates 
of Hq with corresponding eigenvalues Ylrk '’^F\k\nr,Fik) + Y2p^pip)''^pip) — 0 - follows 
that Ho can be identified with a self-adjoint operator on T (cf. Lemma C.3) and that it 
is positive definite. 

Remark 4.3. The significance of the cutoff wq can be explained as follows. The contribu¬ 
tion from the zero mode np(0)^ -|-a;Q^p(0)^ to the free fermion-phonon Hamiltonian is a 
harmonic oscillator Hamiltonian with a discrete spectrum for cuo > 0 but a free-particle 
Hamiltonian with a continuous spectrum for wq = 0. Although the latter is easy to treat, 
with Wo > 0 we avoid a separate discussion of the phonon zero modes which would be 
of no consequence for our results. 0 

^°Some of the notation used here is explained after (2.6). 
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4.3 Bosonized Hamiltonian 


To bosonize the fermion-phonon Hamiltonian in (4.1)-(4.3), we use the boson operators 
arising from the fermion densities and Kronig’s identity; see (2.32), (2.68), and (2.70) 
(note that the latter refers to the free-fermion Hamiltonian Hq in (2.1)). If we use boson 
flavor indices X, X' = F,P and define 


tlpip) = [j+ip) - J-ip)) = (bpiP) - ’ 

^f{p) = . i— (j+{p) + J-{p)) = \f^—r== i^F^P^ + ^U-p)) 

for p ^ 0, then the bosonized Hamiltonian can be written 


(4.11) 




VF {Ql + Q^) + ^g+Q- + : (rip(0)2 + a;24.p(0)2) 


TT 


+ E E Z ^ (fl^(p)n^(p) + u;)^(p)^4>^(p)$^ 

X p^O 

+ E f“np(i5)np(p) + w^(p)^iF(i5)^F 


,0 


p¥=o 

p¥=o 


L 2pvf 


I {p)^p(p) (.p)^Pip) + ^Up)^f 


(4.12) 


where we have collected all zero modes in the first line (recall that Q± = J,j-(0)) and 
introduced 

uj^pip) = vf\p\. (4.13) 

We note that, for both fermions {X = F) and phonons {X = P), the operators ^x{p) 
and tix{p) satisfy the same CCR: 


<l^(p),n^,(p') 


^i^^x,x'5p,pi, 


^Ap),^Up') = ff.(p),rf.,(p') =0 (4.14) 


(this follows from (2.69) and (4.4)), and similarly for the boson creation and annihilation 
operators: 


^xiP)^^^x'iP) =^x,X'Sp,pF bj^{p),b^,{p') 


= 0 , 


b^{p)n = 0. (4.15) 


As a consequence of the results in Sections 2.5 and 4.2, the states 


'/m — 




5t^(0)mp(0) 


RfRZ‘^-n, 


X p^o 


^/mF{p)\ i \/mp(0)! 


mx{p) e No, q± G Z, (4.16) 
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where only finitely many of the mpip) and mp{p) are non-zero, and where m is short 
for the infinite vector {{mF{p))p^ 0 i i'm'P{p))p, Q+, Q-), provide a complete orthonormal 
basis for the fermion-phonon Fock space. These states are also exact eigenstates of the 
free fermion-phonon Hamiltonian Hq (cf. (4.10)) in its bosonized form. 


4.4 Exact solution 

We now turn to the exact solution of the regularized fermion-phonon model. The Hamil¬ 
tonian in (4.12) describes two boson fields with linear coupling terms. Such a Hamilto¬ 
nian can always be diagonalized by a Bogoliubov transformation; see e.g. Appendix C.l 
in [9]. The present case is particularly simple since the computations can be reduced to 
diagonalizing 2x2 matrices (the details can be found in Section 4.6). 

To state our results we find it convenient to introduce the following parameters; 

71 = 72 =-^=, (4.17) 

ZTTVf Vp^TTVf 

which are dimensionless, and 


W = ^J (4 (1 - 72) -vlf + (1 - 71), (4.18) 

which has the dimension of velocity squared. We note that the conditions in (1.4) are 
equivalent to 71 < 1 and 7I < I-I-71, which guarantees that IT, and all other parameters 
given below, are well-defined. We also introduce the operators 


if r = -I-, 




’^r,xiP) — 

This notation is motivated by the following: 


JI,x{p) ^ Jpxi-P) (P>0). 

Jr,F{p) = Jr{p) (cf. (2.32)), and 


(4.19) 


Jr,x{'f’p)^ = 0 Vp > 0, 

Jr,x{p), Jr',X'{p) = r5r,r'Sx,X'-^Sp-p'- 


(4.20) 

(4.21) 


The mathematical properties of the phonon operators Jj.^p[p) are thus identical with 
those of the fermion densities Jr^pip)- For this reason, they can also be used to construct 
phonon vertex operators xexp(i ^p^Q(27r/L)a(p) Jr^p(—p))x (cf. (3.5) for a(0) = tc = 0), 
and the mathematical results obtained for the fermion vertex operators in Section 3.1 
are equally true for such phonon vertex operators. 

We are now ready to formulate a key result towards the exact solution: 
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Proposition 4.4. 

(a) There exists a unitary operator U on T diagonalizing the Hamiltonian in (4.12); 
lA^HIA = £q -\ —z— (Q+ + +Q—) 4“ ^o^p(0)^p(0) 

Jj 


+ ^^v^{p)\p\h'^{p)b^{p) (4.22) 


X p^O 


with 


{vx-vx)\p\^ 


X 0<|p|<'7r/a 

vx{p) = 


vx if \p\ < vr/a, 
vx otherwise, 


(4.23) 

(4.24) 


where^^ 


VF = 


fvl{l-^l)+vl + W 


Vp = 


(1 - 72 ) +vl-W 


(4.25) 


(b) The unitary operator IT commutes with the Klein factors Rr and the fermion charge 
operators Qr- Moreover, 


Px{p) = 


U^Jr{j))U = Ypx{p)jr,x{p) - (Tx{p)J-r,x{p), (4-26) 

ax if 0 < \p\ < tt/o, 


X 


px if 0 <\p\ <Tr/a, 
dx,F otherwise, 


crxip) = 


with the parameters 

fvF 72 Vp{vf + vf{1-7i)) 


PF = 


aF = 



2^^vl-vl (1-7?) 
FF 72t'p {vF - vf{1 - 7i)) 



2VW^vl-vj. (1-7?) 
which satisfy the identity 


PP = 


ap = — 


0 otherwise, 

W l 2 Vp (vp + vf{1 - 7i)) 


(4.27) 


2\/W y'Wp (l - 7?) - Vp 

fvF 72t’P ('Sp - ■vf( 1 - 7l)) 


(4.28) 


2^Jvl [l-^D-v] 


2 2,2 2 1 
Pp — Cp pp — (7p — 1 . 


(4.29) 


(The proof is given in Section 4.6.) 

The formulas for vx, Px-, and ax in terms of the model parameters constitute the 
main computational result. We note that the notation we use for the parameters in (4.28) 
is inspired by [7]. Moreover, the dispersion relations for the fermion-phonon model are 


UJx{p) = vx{p)\p\, 


(4.30) 


^Note that hx > 0 if and only if 71 < 1 and 7 I < 1 -|- 71 . 
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from which it follows that vp and vp have the physical interpretation of renormalized 
fermion and phonon velocities, respectively. One can show, although it is not manifest, 
that our formulas simplify to 


Vp = Vp 



Vp = Vp, 


PF 


Vp + Vp 

2vp ’ 


ap 


± 




Vp — Vp 

2vp 


pp = ap = 0 


(4.31) 


in the limiting case 72 = 0 where the phonons decouple from the fermions; the upper 
sign for ap in (4.31) corresponds to A > 0 and the lower to A < 0. These can be used to 
check that we recover Johnson’s solution of the massless Thirring model [19, 25] (cf. the 
special case Axiaiber = 0 of Klaiber’s solution as in [7], which corresponds to the solution 
by Johnson [21]). 

Remark 4.5. The limit 72 —)• 0 is tricky due to the possibility of incorrectly identifying 
fermions as phonons and vice versa (cf. (4.25)). One way to avoid this is to restrict the 
coupling constants further, as compared to (1.4), and require that (A/(27r))^ < Vp — Vp, 
which corresponds to 7 ^ < 1 — Vp/vp. It then follows that W —)• (l — 7 ^) — Vp in the 
limit 72 —)• 0, which gives the correct identification as in (4.31). 0 

We now turn to the exact eigenstates and eigenvalues of the regularized fermion- 
phonon Hamiltonian. A simple consequence of Proposition 4.4 is the following; 

Corollary 4.6. The exact eigenstates of the Hamiltonian H in (4.1) are given byUrj^ 
with the states in (4.16), and the corresponding eigenvalues are 


^ {q\ + q^_ + 2piq+q-) + LVQmp{H) + EE vx{p)\p\mx{p)- (4.32) 


Thus H defines a self-adjoint operator with pure point spectrum. 

As a special case of this corollary, TL = U8l is the ground state of H with corresponding 
eigenvalue Tq, i.e. Sq has the physical interpretation as ground state energy. 

Proof of Corollary 4-6. It is clear from (2.42a) and (4.22) that HUrj^ = ■ 

Since the states in (4.16) form a complete orthonormal basis in jF, the same holds true 
for these eigenstates. This implies all other results (cf. Lemma C.3). □ 


To compute fermion correlation functions we use the same notation as in (3.24) and 
the regularized fermion field operators in (3.2); see Section 3.3 for motivation and further 
explanations. We thus define 


(O, (xi, fi; 0+) (xnHn] 0+)ij) 

= lim {n,fifi^{xi,ti-,ei).. .if^^{xN,tN]eN)^), (4.33) 

where fJ is the ground state of the fermion-phonon Hamiltonian H and the time evolution 
is given by H, i.e. V’r(x, f; e) = eyip{iHt)fir{x', e) exp(—iJJt). The following result reduces 
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the computation of any such correlation function to an exercise in normal ordering 
products of vertex operators: 


Corollary 4.7. The fermion correlation functions in (4.33) are equal to 

. . .'0^)^(xAr,tAr;O+)O) 

with 

fj'^{x,t-,e) = ^j^qr^-2iTiq((rx-vpt)Qr+'yiVFtQ-r)/L^ 


(4.34) 


Vl 


2tt 


xexp qr'^—px{v)Jr,x{-p) 


X 






(4.35) 


27r 


JJxexp -qr'^—ax{p)J-r,x{-p) 


X 


p^O 




and 


2tt 


Za,e = exp I - ^ ^ I ■ 

0<p<7r/a 


(4.36) 


(The proof can be found further below.) 

By arguments similar to those in Section 3.3, it follows from the corollary above that 
the correlation functions in (4.33) are equal to 

h 7 \ 

^ n (4.37) 

\\TJ i<n<m<N 


with building blocks C^'^, (•,•;•) which can be computed by normal ordering products of 
two operators of the form in (4.35): 


Vi^(x, t; €)i>l, {x', t'; e') = Cf!^, {x - x ,t - t'; e + e') iip'^ix, t; e)^^, (x', t'; e')> 


(4.38) 


In Section 4.5, we give explicit formulas for these building blocks and for all fermion 
correlation functions in the continuum and thermodynamic limits a —)• 0“*" and L —)■ oo. 
Fermion two-point correlation functions for a > 0 and L < oo can be found in [28]. 

Proof of Corollary f.l. To obtain (4.34) we insert 0. = UCl and I = 1414^^ into (4.33) and 
introduce 

■0^(x, t] e) = I4^if^{x, t; e)l4 = e*^*Vj'?(x; e)e“*^* (4.39) 

with H = hf^HU in (4.22). The explicit form of the operators in (4.39) can be computed 
using (3.20)-(3.21) and Proposition 4.4 as follows. Since 14 commutes with Rr and Qr-, 


fj^ixie) = 


K 

\/2xe 


2tt. 


^i?re-2"^"^<5-/^xexp I | , (4.40) 
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where it follows from (4.26) and (4.21) that the rightmost exponential is equal to 


2tx 


exp Qr^—px (p) Jr,x (-p) 


X 


p^O 




27r 


exp -qr ^ —ax{p)J-r,x{-p) 


X 






(4.41) 


We are left to compute the time evolution given by H: for the density operators, 

e^^^Jr,x{p)e~"^' = ^r,x(p)e-™(^’)*, (4.42) 

and for the Klein factors, 

= ^R^e^i^rvp(Qr+'nQ-r)t/Lx _ 

Combining these with (4.39)-(4.41) yields the result in (4.35); all that remains is to write 
this expression in a boson-normal-ordered form using (3.11a) and (3.9), which gives rise 
to the factor 


\/1 — gxp 

(this follows from (3.21), (4.29), and (4.36)). 


(Pf + O-F + /?P + CTp - 1) 

0<p<7r/a 


-ep _ 


Me 


Za 


(4.44) 

□ 


4.5 Continuum limit 

As mentioned after Corollary 4.7, all fermion correlation functions in the fermion-phonon 
model can be computed from the results in Section 4.4. These formulas are in general 
complicated, but if one is interested in intermediate length scales they can be simplified. 
One way to compute such simplified formulas is as limits a —)■ O'*' and L —oo of the 
fermion correlation functions in the regularized model; this is, however, demanding from 
a computational point of view. In this section we construct a renormalized model which 
corresponds to the continuum limit a —)■ O'*" of the regularized fermion-phonon model. 
Using this renormalized model, the continuum-limit results for the fermion correlation 
functions can be obtained with less computational effort, but, as will be discussed, the 
continuum limit requires additive and multiplicative renormalizations. 

To see that the limit a —)■ 0“*" is non-trivial, note that the ground state energy in 
(4.23) diverges like OiLjc?^. Moreover, for e = 0, the constant in (4.36) is vanishing in 

^^This can be shown by differentiating the l.h.s. of (4.42) with respect to t, computing the commutator 
of Jr(p) with H =U^HLJ by recalling (4.19), and by integrating the resulting differential equation. 

^^This can be computed using (2.38) for p = 0; the computation is similar to that of the time evolution 
of the density operators. 
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this limit: 


Za,o — 





(l + O 



(4.45) 


with the Euler-Mascheroni constant 7 = 0.5572.. The reason for these difficulties is 
that the unitary operator U in Proposition 4.4 does not exist in this limit. It is therefore 
remarkable that all these divergencies can be eliminated in a simple way after a similarity 
transformation with U using additive and multiplicative renormalizations as follows. We 
define a renormalized Hamiltonian 


Hj-en = lim 
a—>-0+ 




-£c 


(4.46) 


and renormalized fields 

^^{x, t; e) = lim I — ) t; e)U, (4.47) 

a^ 0 + ' \ L J 

where £ > 0 is an arbitrary length parameter introduced for dimensional reasons.It 
follows from the results in Section 4.4 that these define a renormalized fermion-phonon 
model in the following sense: 

Proposition 4.8. 

(a) The Hamiltonian in (4.46) and the fields in (4.47) are equal to 

Hren= {Q+2-/iQ+Q-)uob^p{0)bp{0)ujx {p)bxip)bx(.p) ( 4 - 48 ) 


and 


T^(x,f;e) = 


1 /27ri 

7 t 


^^j^qr^-2'Kiq({rx—vpt)Qr+'iivptQ-r)/L>^ 


X n^exp I qrY^‘^p^X^^{-p)e-^P^-''-xt)-e\p\/2 

X \ p^O 


X TT xexp I -qrS^ ‘^axJ-pxi-p) 

X I U 




(4.49) 


and they define self-adjoint and unitary operators on T, respectively. 

(b) The renormalized fermion correlation functions 

f n\N{(j\+a-y) _ 

lim lim ( e'^—) (11, V'®) (xi, ti; 0+)... V’^)^(xAr, fAr; 0+)n) 

L^>ooa-:>0+ \ a J 

^"^This follows from 1/u = 7 + In(Af) + Ofl/X). 

^®Note that Z~f Vf’’ —>■ (l + O (a/L)) in the limit e —>■ 0^. 


(4.50) 
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exist and are identieal with 


lim (4.51) 

L^oo 

as funetions away from singularities}^ 

(The proof can be found at the end of this section.) 

Remark 4.9. We only prove that the correlation functions in (4.50) and (4.51) are iden¬ 
tical as functions away from the singularities. To see why, we emphasize that the second 
expression is shown to equal 

/ f\N{al+aj,) _ 

lim lim lim i e^ — ) (fl, (xi, ti; ei)...'0^);)(xAr, tAr; eAr)0) (4.52) 

L —^CXD —^0"^ CL —^0*^ y ^ / 

in the proof of Proposition 4.8, and this differs from the first by the order in which the 
limits ej —)■ O"'' and a —)■ 0+ are taken. We cannot exclude the possibility that these 
limits cannot be interchanged, and, if that were the case, then the correlation functions 
in (4.50) and (4.51) would define different distributions.^^ We emphasize that, from a 
physical point of view, the correct order of limits is as in (4.50). 0 

As for the regularized model, the renormalized Hamiltonian P^ren determines the time 
evolution of the renormalized fields, i.e. 'kr(x,t;e) = exp(zijreni)'kr(x; e) exp(—iiJren^)- 
It would be interesting to extend this renormalized fermion-phonon model to allow for 
an efficient computation of boson correlation functions. Another interesting use would 
be to compute correlation functions at finite temperature; see [9] for such results for a 
similar model. However, this is beyond the scope of the present paper. 

Let us briefly discuss the physical interpretation of the renormalized fermion-phonon 
model. Our construction of the regularized fermion-phonon model is on the Fock space 
for the corresponding non-interacting model, and the Fock vacuum H has the physical 
interpretation as ground state of the latter model. The ground state H of the interacting 
model with local interactions (which corresponds to the limit a —)■ O'*') does not exist 
in this Fock space; this is clear since, for example, the ground state energy diverges. 
However, using a similarity transformation with the unitary operator li, one can con¬ 
struct a new Hilbert space which is such that, by definition, it contains the model ground 
state: the definition is such that the transformation maps Q. =hlCl to the Fock vacuum 
H Since the normal ordering was defined with reference to the original Fock vac¬ 

uum, the additive and multiplicative renormalizations can be interpreted as changing 
the normal ordering such that it uses the interacting ground state as reference instead. 
The difference between these two normal orderings diverges in the limit a —)■ O"*", but 
this is to be expected since the similarity transformation maps the original Fock vacuum 
to the state which does not exist in this limit. 

^®This is the strongest statement of this identity which we prove in this paper; see Remark 4.9. 

^^This is analogons to the situation for l/{x± iO"*") and "Pfl/x) which are identical as fnnctions away 
from the singularity in a; = 0 but define different distributions. 
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The renormalized model can be used to efficiently compute the continuum and ther¬ 
modynamic limits a —7- 0^ and L —)■ oo of the renormalized fermion correlation functions. 
The strategy is the same as described for the regularized model at the end of Section 4.4. 
The computational advantage of working with the renormalized model is that each of the 
renormalized fields Tr(x,t;e) in (4.49) is essentially a product of four vertex operators 
of the form 


e) = xexp v 


27r - 


-p)e 


—ipx—e\p\/2 


(4.53) 


with different real parameters v (cf. (4.49)). We also have the following general formula 
for writing products of such vertex operators in a normal-ordered form: 


W,^^(x;e)W,^:^( 


/ i exp (- f[r{x- x') + 5 (e + e')]) \ 
2sin(f[r(x-x0 + |(e + e')]) ) 


(4.54) 


with X- • -x meaning boson normal ordering of the l.h.s. and 

2ti(. iexp (-^[rx-k |e]) i(. 

L 2sin (j^[rx-|-jc]) rx-|-f0+ 


(4.55) 


(this follows from (3.14a) using (A.l)). Moreover, it is useful to note that the zero-mode 
contributions when normal ordering products of the renormalized fields Tr(x, t; e) can be 
ignored in the limit L —)■ 00 (since they always come with a factor 1/T); see Remark 4.11 
for a discussion of the zero-mode contributions for L < 00 . It thus follows that 


T^(x, t; e)TS(x',t'; e') = (x - x',t- t'; e + e) x'i;'^(x,t; e)'I'S(x', t'; e')x (4.56) 


with 


/27r.^\ 


C±,±{x,t;e) 


e-:>0+ -*■ i \rx — Vxt + f0+ J 

r,X 


- ' C’±,^(x,t;e) —^ II 


L J 


U 


e^0+ V ~ Vxt + 

r,A 


\ px(yx 


(4.57) 


As we mentioned after Corollary 4.7, by arguments similar to those in Section 3.3, this 
implies the following: 


Corollary 4.10. In the thermodynamic limit L —)■ 00 , the renormalized fermion corre¬ 
lation functions in (4.51) are 


lim (fl,T^)(xi,ti;0+)...^^));(xAr,tiv;0+)fI) = {H, 


L —^CXD 


1 


^ .27r£ ) 


n n 

0<n<m<N r,X 


a 


r{Xn - Xm) - Vx{tn “ ^m) + ^0+ 


QnQmCr^X;rn,rm 


(4.58) 
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with 


Cr,X-,±,± — Px^r,± + Cr,X;±,=F = PXO'X, 


where the first factor on the r.h.s. of (4.58) is a sign factor {cf. (3.26)ff). 


(4.59) 


Proof. The formulas for Cr^X]rn,rm follow from the exponents in (4.57), and we are thus 
left to prove that the factors (27r£/L)°'f’+°'p from the renormalized fields in (4.49) exactly 
match the factors needed for the limit L —)■ oo to exist. This is equivalent to showing 
that 


1 f27re\ 


0-p-|-(Tp\ ^ 


/j_xiV /2 

V27r.^y 


n 

l<n<m<N 


V ^ J 

if (n, 0, which is true by Corollary 3.11. 


P 0'TrP\ 

‘ 1 (4.60) 


□ 


The explicit expressions for all fermion correlation functions can be found from Corol¬ 
lary 4.10. For example, the special case N = 2 gives the two-point correlation functions 
stated in the introduction (cf. (1-5)). Similarly, as an example of a four-point correlation 
function (N = 4), we find 


lim (n, T+(xi, ti; 0+)4 'L(x 2, t 2 ; 0 +)T_(x 3 , ts; 0 +)TY(x 4, t4; 0+)^) 

L—>oo 

_ / J_\ ^ "rr ^ rxi3 - 7x^13 + / rX24 - + i0''~ \ 

\2tt£ J ^}l\rxi 2 - Vxtl 2 + iO'^ J V™34 -'Sxi34 + iO+y 


n 

X 

n 

X 


a 


3^14 — VxtlA + i0+ 
i£ 

X23 - Vxt23 + i0+ 


Px 


i£ 


Px^x 


A 


—Xi4 — 7x^14 + *0+ 

_ u _ 

-X23 - Vxt23 + , 




using the shorthand notation Xnm ^ Xn — Xm and tnm ^tn — tm- 

Remark 4.11. One can get simple formulas even without passing to the thermodynamic 
limit: for the fermion two-point correlation functions we checked that the zero-mode 
contributions exactly cancel the (x ± ri;xi)-dependence from the numerator in (4.54) by 
using the identity 

{p‘f + xf + {p'p + PPp) vp = vp, (4.62) 

which can be derived from (4.28). We expect that a similar cancelation takes place for 
all fermion correlation functions. 0 


Proof of Proposition 4-8. The states in (4.16) form a complete set of orthonormal eigen¬ 
states of the renormalized Hamiltonian in (4.48), and the corresponding eigenvalues are 
lim„_,.g+(Tm — To) with Tm in (4.32). It follows that Hren defines a self-adjoint operator 
on R (cf. Lemma C.3). This also shows that the limit in (4.46) is in the strong sense on 
the domain of finite linear combinations of the states in (4.16). 
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4 Fermion-phonon model 


Equation (4.47) implies (4.49) for the renormalized fields. They are thus proportional 
to a product of unitary operators (this follows from Proposition 3.4) and are therefore 
also unitary. This completes the proof of Part (a). 

The limit in (4.47) is in the weak sense on T", and it thus follows that the correlation 
functions in (4.51) and (4.52) are equal. This proves Part (b) (see Remark 4.9 for further 
explanations). □ 


4.6 Proof of Proposition 4.4 


One important feature allowing to solve the fermion-phonon model exactly is that bosons 
with different magnitude of momentum p decouple. To make this manifest we write the 
Hamiltonian in (4.12) as 

H = '^hp (4.63) 

p>o 

(this follows since {hp + h-p)/2 = hp) with the zero-mode part 

/iQ = ^ (vp {Q+ + Q- + “^nQ+Q-) + np(0)^ -l- a;Q4>p(0)^^ — -ujo (4.64) 

and the parts with momenta p > 0 written as 

/ip = ^ (n^(p)A(p)n(p) ^\p)B{p)^{p)'^ - {u^pip) + u}p{p)) (4.65) 

using the following matrix notation; 


h{p) = {Up{p),Up{p)) 

A(.). ;) 


T 


$(p) = (l>p(p),4>p(p))^, 


VpVp-f2{p)\ 

4 J 


(4.66) 

(4.67) 


with the shorthand notation 7i,2(p) = 7i,20(7r/a — |p|). To ensure that this describes a 
stable system, we must require that both matrices A(p) and B(p) are positive definite. 
This is true if and only if 71 < 1 and 7I < I-I-71, which, as we mentioned, are equivalent 
to the conditions in (1.4). 

The parts with p = 0 and p > pja are already diagonal: 


ho — -^—{Q\ + Q^_ + 2^iQpQ-) + a;o6p(0)6p(0), (4.68) 

hp>-n/a = + b^xi-p)^x(-P)) ■ ( 4 - 69 ) 

X 


The task of diagonalizing the parts /io<p<7r/a can be reduced to diagonalizing the matrix 

C(rt . A,p)./^B(p)A(p)V. = 7 , p™. 
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which is also positive definite if and only if 71 < 1 and 72 < 1 + 71 • It follows that there 
exists a unitary matrix U and a diagonal matrix u{p) = diag{u f{p), up(p)) such that 

C{p) = Uu{pf\jK (4.71) 

Indeed, one can write 

ho<p<^/a = Y (n^p)n(p) ^\p)u{pf^{p)'^ - {u%{p) + U}p{p)) (4.72) 

with 

n{p) = \JA^/^U{p), $(p) = U'^A-^/2 ^(p) (4.73) 

defined in the same way as in (4.66). This can be shown to be a canonical transformation, 
i.e. one can show that the operators Ilx{p) and ^x{p) satisfy CCR similar to (4.14). By 
introducing boson operators 6^^(p) corresponding to the ones in (4.9) and (4.11); 

^xiP) = " ^Xi-P)) > 

I — (4-74) 

one obtains 

ho<p<nla = Y1 {^x(p) (b^xipy^xip) + '^xi-pyxi-p)) + {‘^x{p) “ Wx(p))) • (4-75) 

X 

It is known that there exist unitary operators Up such that 

Ul,by\±p)Up = by\±p), (4.76) 

which can be constructed using the operators 6^^(±p); see e.g. Appendix C.l in [9] for 
an explicit construction of these unitary operators. It follows that the unitary operator 

U= II Up (4.77) 

0<p<7r/a 

satisfies (4.22)~(4.23) (note that this is a finite product and thus well-defined). Moreover, 
one can easily verify (4.24)-(4.25) by computing the eigenvalues of the matrix in (4.70) 
(cf. (4.30) and (4.71)). This completes the proof of Part (a) of Proposition 4.4. 

To prove Part (b), invert (4.73) and write the resulting equations in component form; 

IIx(p) = ^ i^n)x,x' ^x'{p), 4*x(p) = ^ {^^)x,x' ^x'{p) (4.78) 

X' X' 

using matrices Mn = and M$ = A^/^U (as before, we restrict our discussion 

to the non-trivial case 0 < |p| < p/a). Using (4.9), (4.11), and (4.74) one finds 
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with 


Cx,x' - 2 ( \[^,{^^)x,X' ^\l i^^)x,x' ) ’ 


VX' 

vx 

VX' 

vx 


= 2 l^V■ V ) ■ 

It thus follows from (4.76)-(4.77) that 

^^b^{p)U = (^^x,x'bx'(p') '^x,x'b\c'(~P)J ’ 


and recalling the relation between b^x\p) Jr,x{p) in (4.19) we obtain 
Jr,X{P)U = E( Cx,X'Jr,X'{p) + Sx,X'J-r,X' 


(4.80) 


(4.81) 


(4.82) 


which, for X = F, implies (4.26)“-(4.27) with px' = Cf,x' and ax' = —Sf,x'- The 
explicit expressions in (4.28) follow from (4.80) by straightforward computations which 
we omit. Lastly, the consistency of (4.81) with the CCR in (4.15) implies 


(4.83) 


= 1 


X' 


which, for X = F, is equivalent to the identity in (4.29). 


5 Concluding remarks 

We constructed and solved a one-dimensional fermion-phonon model using bosonization. 
We also gave a self-contained account of the mathematical basis of bosonization together 
with complete proofs. Our main results include a renormalized model in the continuum 
limit, which was obtained by additive and multiplicative renormalizations (cf. Proposi¬ 
tion 4.8), and explicit formulas for all fermion correlation functions in the continuum and 
thermodynamic limits using this renormalized model (cf. Corollary 4.10). We emphasize 
that positivity of our correlation functions is manifest from the method by which they 
were computed: the renormalized model we derived is analogous to the one used by 
Carey, Ruijsenaars, and Wright [7] to prove positivity of Klaiber’s X-point correlation 
functions for the massless Thirring model [21], and their argument also applies to our 
model. 

As an application of the tools developed in this paper, we address one point which 
was raised in the introduction: we mentioned that certain fermion four-point correlation 
functions can be used to study physical properties of the fermion-phonon model related 
to charge-density-wave (CDW) and superconducting (SC) instabilities [8, 40]. The corre¬ 
lation functions of interest are lim2,_^oo(f^, Oyl(2:T)OA(0) 0)11) for A G {CDW, SC} with 
OcT)w{xit) = x'I'_,_(x, t; 0''')'k[L(a::, t; 0''')x andOsc(a:,f) = xT_ (x, t; 0''')T_,_(a:, t; 0+)x . 
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We emphasize that these operators must be in a boson-normal-ordered form for the cor¬ 
relation functions to exist, and this corresponds to an additional multiplicative renormal¬ 
ization; this is seen by comparing the equations below with (4.61). Using Corollary 4.10 
we find 


lim (n, O 

L^oo 


CDW 


(x, 


( 0 , 0 ) 0 ) 



e 

(vxt — t0+)2 


{px-^xY 


(5.1) 


and 


lim (0, Ogc (x,t)OgQ(o, o)r2) 

L—>-oo 



e 

[vxt — i0+)2 




(5.2) 


which are of interest in physics since they allow to construct phase diagrams [8, 40]. 

We derived our results using a particular regularization, but, as we mentioned, there 
exist a vast number of different regularizations. Physical arguments suggest that the 
correlation functions in the continuum and thermodynamic limits are universal in the 
sense that they are independent of the details of the regularization. It would be in¬ 
teresting to prove this universality conjecture. A related question is if the correlation 
functions obtained by passing to the continuum limit are identical as distributions to the 
continuum limit of the correlation functions for the regularized model. (We only proved 
the equality of these functions in a weaker sense.) As for the previous question, physical 
arguments suggest that this is the case, but to prove this conjecture one would have to 
show that the correlation functions in (4.50) and (4.51) are identical as distributions. To 
our knowledge, both of these questions are open even for the massless Thirring model. 


Acknowledgments: We would like to thank Alan Carey, Jan Derezihski, Luca Fresta, 
and Jonas de Woul for valuable suggestions on the manuscript, and Vieri Mastropietro 
for encouraging us to write the paper. Financial support by the Goran Gustafsson 
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Appendix A Notation and conventions 

The symbols C, M, and Z denote the sets of complex, real, and integer numbers, and 
N and No denote the positive and non-negative integers, respectively. The symbol dx is 
short for djdx. We write “=” to emphasize a definition. The complex conjugate of c G C 
is c. The Heaviside function ©(•) is defined such that 0(x) = 1 for x > 0 and ©(x) = 0 
for X < 0. We use the symbols [a, b] = ab — ba and {a, b} = ab + ba for the commutator 
and the anticommutator, respectively. We write I for the identity operator, but often 
identify c G C with cl. Two symbols f and * are used for Hermitian conjugation of 
operators: the first is used in cases where operator domains are well-defined or of no 
concern, and the second in a weaker sense such that, if A is an operator defined on a 
dense domain "D, then “Al” means “A* restricted to U.” 
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Appendix A Notation and conventions 


Unless otherwise stated, and when there is no risk of confusion, a free variable or 
index in an identity can take any allowed value. For example, the relations in (2.2) hold 
true for all k,k' G (27r/L)(Z +1/2) and r, r' = ±, whereas (2.3) holds true for r = ± but 
only those k G (27r/L)(Z + 1/2) which satisfy the given condition. 


Parameters: The following parameters are used: 


Parameter 

Physical significance 

L 

infrared cutoff (length of space) 

a 

ultraviolet cutoff (interaction range) 

VF 

Fermi velocity 

Vp 

phonon velocity 

A 

fermion-fermion coupling constant 

g 

fermion-phonon coupling constant 


Variables: The following variable names are reserved for particular sets; 


Variables 

Set 

Physical interpretation 

r, r', rj 

{+) “} 

chiralities 

x,x',y,xj 

[-L/2,L/2] 

positions 

k, k\ kj 

(27r/L)(Z + l/2) 

fermion momenta 

p,p',Pj 

{2ttIL)'L 

boson momenta 

X, X' 

{F,P] 

boson flavor indices 


Summation conventions: The following abbreviations are used: 


Abbreviation 

Meaning 

E 

E 

r 

r=± 

E 

E 

k 


E 

E 

p¥=o 

pe^z\{o} 

(etc. 

) 


Useful formula: The following formula is used repeatedly: 

= (Re(^) > 0). (A.l) 

p>0 
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Appendix B Bosonization in a nutshell 

In this appendix we give a survey of the results developed in this paper using a notation 
common in the physics literature. It serves to connect the results in this paper to the 
latter and as a summary. We also include remarks on the history of bosonization. 

B.l Bosonization 

The collection of mathematical results known as bosonization concern the quantum held 
theory model dehned by the Hamiltonian 

/.L/2 

Ho= dxvp'^ ■■ipl{x)r{-idx)'ipr{x)- (B.l) 

J-L/2 j. 

with fermion helds satisfying antiperiodic boundary conditions and the CAR 

= ^r,r'S{x - y), {(x), (y)} = 0. (B.2) 

The fermion densities are dehned as 

Jr{x) = :^l{x)il!^{x): . (B.3) 


By bosonization we mean the following three identities: 


The fermion densities obey non-trivial commutation relations: 

[Jrix], Jr'iv)] = r^drydxS{x - y). 

Zm 


(B.4) 


The Hamiltonian can be expressed in terms of fermion densities: 

rL/2 


j-n/z 

Hq = 'kvf\^ / dx ■.Jr{x)‘^ . 
„ J-L/2 


(B.5) 


The fermion fields can be expressed in terms of fermion densities: 


iprix) = exp {r2'Kid^ ^ Jr (x)) . 


(B.6) 


The r.h.s. of (B.4) is an example of a quantum held theory anomaly known as a Schwinger 
term, and the result in (B.5) is known by the name Kronig’s identity}^ Moreover, the 

Actually, what is usually called Kronig’s identity is 



da; :‘ipl{x)r{—idx)ipr{x)- 


xL/2 

J-L/2 


d® :Jr{x): ■ 


However, we only need the somewhat weaker form stated in (B.5). 
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Appendix B Bosonization in a nutshell 


inverse differential operator d~^ in (B.6) hides an important integration constant, which 
includes an infinite normalization constant and a so-called Klein factor (cf. (3.20)). Note 
that all three identities are consequences of mathematical subtleties related to normal 
ordering. 

To explain the name bosonization we note that 

nir(x) = -^UF7r[j+(x) - J-(x)], $F(a:) = 5“^/^[j+(x) + J_(x)] (B.7) 

V 

are boson field operators since they are Hermitian and satisfy CCR: 

[^F{x),IlF{y)] =iS{x-y), [^>f(x), ^^(y)] = [ni;’(x), nF(y)] = 0, (B.8) 

and it follows from (B.5) that the fermion Hamiltonian in (B.l) is identical with a boson 
Hamiltonian: 

= 2 y dx : (ni.(x)2 + 4 [aa.^>F(x)]2 j :. (B.9) 

Finally, the identity in (B.6) can be used to express all fermion correlation functions in 
terms of boson correlation functions, and the latter can be computed exactly. 

Remark B.l. By interpreting vf as the vacuum speed of light, the Hamiltonians in (B.l) 
and (B.9) define important models in particle physics: the former corresponds to massless 
Dirac fermions, and the latter to massless Klein-Gordon bosons, both in 1+1 dimensions. 
The inverse differential operator in (B.7), however, hides certain complications related 
to zero modes; see the closing discussion in Section 2 for details. 0 

Remark B.2. From the point of view of particle physics described in the previous remark, 
the fermion densities in (B.3) correspond to fermion currents in light-cone coordinates. 
This explains why we use the symbol Jr{x) for these operators. 0 

B.2 Formal solution of the fermion-phonon model 

The formal bosonization results in (B.1)-(B.9) can be used to bosonize the fermion- 
phonon model in (1.1)-(1.4) as follows. By restating the model with an infrared cutoff 
L > 0, we note that the free-fermion Hamiltonian in (B.l) is part of the fermion-phonon 
Hamiltonian with the cutoff in place. As a result, this entire Hamiltonian can be cast in 
the form 

rL/2 /i 1 

dx I -:(ni7’(x)^+ 4[5a:^F(a:)]^): +-:(np(x)^+4[5x^p(a;)]^): 

J — Z//2 \ 

+ {-IlF{xf + VF[da;^F{x)f) + g^^^[dx<^F{x)][dx^p{x)^ (B.IO) 

consisting only of bosons field operators 4>x(a;) and nx(x) in the sense that the latter 
are Hermitian and satisfy CCR: 

^x{x),Ilx'iy) =i5{x-y), ^x{x),^x'iy) = Iix{x),Ilx'{y) =0. (B.ll) 
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A Hamiltonian of this type can be diagonalized by a Bogoliubov transformation. More¬ 
over, as in the free-fermion case, all fermion correlation functions can be computed using 
the identity in (B.6). In the formal discussion given here, we have ignored serious math¬ 
ematical difficulties. For instance, it is standard practice to construct a Hilbert space for 
the non-interacting model, i.e. the Hamiltonian in (1.1) with X = g = 0. However, this 
Hilbert space does not contain the ground state of the interacting model. Furthermore, 
one can show that the exact fermion two-point correlation functions for the interacting 
model are inconsistent with the CAR in (B.2). This was discussed by Wilson in [42] 
for the massless Thirring model, which corresponds to the special case <7 = 0, and it is 
equally true for gr / 0. 

B.3 From formal to precise results 

To avoid difficulties with the distributional nature of quantum fields, we work with field 
operators in Fourier space: 

1 

Tprik) = / dx'ip^{x)e~^ Jr{p) = / dx J,., (B-12) 

V27r J-L/2 J-L/2 

where k G (27r/L)(Z -|- 1/2) and p G (27r/L)Z are fermion and boson momenta, respec¬ 
tively. (Note that one can regard and Jr{p) as fermion field operators and density 

operators, respectively, smeared out by test functions, the latter being plane waves.) This 
can be used to write down and prove precise versions of the formal bosonization results 
above; see Sections 2 and 3. Similarly, as in the case of free fermions, our strategy to 
make the fermion-phonon model precise is also to work in Fourier space: 

Up{p) = -= dxUp{x)e-^P^, 4>p(p) = ^/ dx<Pp{x)e-^P^, (B.13) 
V27r J-L/2 J-L/2 

1-1/2 i-L/2 

^a{p)=l dxAa(x)e"*^’*, g^{p)= dx ga{x)e~"P'^. (B.14) 

J-L/2 J-L/2 

This and (B.12) can be used to obtain the precise formulation of the fermion-phonon 
model in Section 4 from the formal definition in the introduction. 


Appendix C Operators and quadratic forms 

Many non-trivial mathematical questions in quantum physics are due to unbounded 
(Hilbert space) operators. For instance, products and commutators of unbounded oper¬ 
ators need not be defined, formally symmetric operators need not be self-adjoint, etc.; 
see e.g. [30]. To keep this paper self-contained, we collect a number of results which 
allow us to settle such questions for the operators encountered in the main text. 
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C.l Quadratic forms 

Consider a complex Hilbert space T with inner product (•,•), and let P be a dense subset 
of T. A quadratic form g on P is a map "D x "D —)■ C, (/, ff) e-)■ q(f, g) such that q{-,g) is 
conjugate linear and q{f, •) is linear; see e.g. Section VIII.6 in [30]. 

Choose some fixed complete orthonormal basis in iF labeled by a countable set of 
quantum numbers; we write this basis as without loss of generality. We are 

mainly interested in the case where is a fermion Fock space and {j^nlneN either 
the fermion states in (2.6) or the boson states in (2.56). 

Let V be the set of all finite linear combinations of the basis states rjn,^^ and let 


A — {An,m)n,m£N (C-1) 

be some infinite matrix with complex entries. Such a matrix can always be identified 
with a well-defined quadratic form on T> as follows; 

OO 

qA{f,g)= X] An,m{f,rin){rim,g) {f , d ^ T^) (C.2) 

n,m=l 

(note that this sum is always finite). Moreover, the adjoint of A, given by 


(A^)m,n 


= A. 


n,m^ 


(C.3) 


is always well-defined as a quadratic form. However, to be able to identify this quadratic 
form with an operator on iF: 


OO 

Af ^ ^ An,rar]n{Vm, f) (C.4) 

n,m=l 

one has to impose certain restrictions on the matrix elements An^m- Note that we abuse 
notation and use the same symbol A for the matrix in (C.l), the corresponding quadratic 
form in (C.2), and, if it exists, the corresponding operator on F'. If a quadratic form A 
can be identified with an operator, then the matrix elements correspond to expectation 
values of this operator: 

An,m — {jlm Ag^n) • (^. 5 ) 

It should be noted that the product of two quadratic forms A and B need not be a well- 
defined operator even if both forms can be identified with operators. However, under 
certain restrictions on the matrix elements, it is true. In such cases, the operator product 
corresponds to the matrix product; 

OO 

An,kBk,m- (^.6) 

fc=l 


^®To be precise, T) is the set of all states g — CnVn with c„ £ C and finite N G N. 
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C.2 Four lemmas 

We now define a set of quadratic forms as in (C.2) satisfying conditions which rule out 
many of the possible difficulties mentioned. 

Definition C.l. Let Of be the set of all matrices A as in (C.l) such that there exists 
an M < oo such that, for all n G N, 

An,ra = Ara,n = 0 \fm> M. (C.7) 

From this definition it follows that A £ Of if and only if Atj, A^r] £ V for all rj £1). 
The following result is therefore obvious: 

Lemma C.2. The set Of can be identified with a *-algebra of operators on T with 
involution f and identity I {i.e. In,m = dn,m)- 

The above lemma makes precise in which sense many of the products and commu¬ 
tators of operators in the main text are well-defined. 

It should be noted that an operator A £ Of satisfying A = A^ does not necessarily 
define a self-adjoint operator; see e.g. [16] for the stronger conditions needed to ensure 
this. However, there is one important special case where this is trivially true: 

Lemma C.3. Given A £ Of such that 

Ar]n = (inT]n with On G M Vn G N. (C.8) 

Then A is the restriction of a unique self-adjoint operator on if. 

(This is a simple consequence of the definition of self-adjoint operators; see e.g. [30].) 

We often abuse notation and use the same symbol for the self-adjoint extension of 
such an operator A. 

In the main text we introduce vertex operators as quadratic forms of the kind dis¬ 
cussed above, and we show that the product in (C.6) of two such quadratic forms is well- 
defined. There is an important case when one can conclude from this that a quadratic 
form can be identified with a unitary operator: 

Lemma C.4. Let A be a quadratic form as defined in (C.1)-(C.2) and such that A^A 
and AA^ both are well-defined and equal to the identity I. Then A can be identified with 
a unique unitary operator on F. 

If this is the case we abuse notation and use the same symbol for the unitary operator 
extending A. 

Proof. Let f £ F and f^^\ n G N, a Cauchy sequence in F converging to /. Then A'^A = 
I implies that Af^'^'^ is also a Cauchy sequence in V with a unique limit lim„_^oo Af^^'i = 
Af depending only on /. In a similar manner, AA^ = I implies that A^ has a unique 
extension from F to F. It is easy to see that these extensions satisfy A^ A f = AA^^ f = f 
for all f £F. □ 
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In the main text we use a simple strategy to establish commutation relations of 
unbounded operators with a common dense invariant domain V>: we introduce a cutoff 
A to approximate the unbounded operators by bounded ones (for which domain questions 
can be ignored), and we derive the result by taking the limit where the cutoff is removed; 
see also Remark 2.4. For the convenience of the reader we state a definition and a lemma 
underlying this strategy. 

Definition C.5. (Operator convergence) Suppose A ^ Of and let {yl^}A>o be a 
family of bounded operators. Then 

A = lim A^ on V (C.9) 

A^OO 


means the following: for fixed / G P, the identity A^f = Af holds independent of A for 
sufficiently large values (i.e. 3Ao < oo such that A^f = Af VA > Aq). 

Lemma C.6. Let {A^}a>o be a family of bounded operators and suppose A G Of such 
that (C.9) holds true, and similarly for {R"^}a>o 0 ,'i^d B {if B is bounded one can set 
B^ = B). For fixed f gT>, assume there exists some C G Of such that 




f 


Cf 


(C.IO) 


holds true independent of A and A' for sufficiently large values. Then [A, B] = C onV. 
(The proof is straightforward and thus omitted.) 
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